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B.l  Goodness  of  fit  test  data  for  ( i,j ) = (10,22)  member  of  the  parameter- 
ization set.  Column  e(i)  denotes  expected  values,  while  y(i)  denotes 
measured  values.  The  chosen  parameterization  was  A#  = 7r / 64  and 
Ap  = h 135 

B.2  Goodness  of  fit  test  data  for  (i,j)  — (55,  7)  member  of  the  parameter- 
ization set.  Column  e(i)  denotes  expected  values,  while  y(i)  denotes 
measured  values.  The  chosen  parameterization  was  AO  = 7t/64  and 
Ap  = h 136 

B.3  Goodness  of  fit  test  data  for  (i,j)  = (200,22)  member  of  the  parame- 
terization set.  Column  e(i)  denotes  expected  values,  while  y(i)  denotes 
measured  values.  The  parameterization  was  AO  = 7r/64  and  Ap  = h.  137 
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The  most  popular  method  for  the  recognition  of  straight  lines  in  a digital  image  is 
the  Hough  transform.  In  this  work  we  study  the  effects  of  the  digitization  errors  on  the 
Hough  transform  with  the  tools  provided  by  the  Digital  Topology  Theory.  We  derive 
a corrected  formula  for  the  worst  case  spreading  of  a digital  straight  line  segment  of 
a given  length.  We  also  show  that  even  in  the  perfect  case  (no  noise,  perfect  digital 
straight  line  segments)  the  peak  cell  may  not  best  represent  the  parameters  of  the 
digital  straight  line  segment.  We  also  study  the  performance  prediction  problem  and 
prove  that  if  the  noise  in  the  image  is  uniform  the  number  of  votes  in  the  parameter 
space  follows  a binomial  distribution.  We  use  the  theory  presented  to  propose  an 
algorithm  for  the  estimation  of  the  length  of  a digital  straight  line  segment  from 
its  Hough  transform.  The  algorithm  is  extended  to  cases  where  noise  and  other 

objects  exist.  Based  on  the  length  estimation  algorithm,  we  propose  a new  algorithm 
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for  the  detection  of  convex  polygons  in  a digital  image.  The  method  is  tested  in 
both  synthetic  and  real  images  and  gives  very  satisfactory  results.  We  also  prove 
the  uniqueness  of  the  representation  of  a convex  polygon  from  the  peaks  it  gives 
in  the  Hough  (parameter)  space.  Finally,  we  propose  an  efficient  algorithm  for  the 
construction  of  a convex  polygon  from  the  parameters  of  its  edges. 
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CHAPTER  1 
INTRODUCTION 

1.1  Purpose  of  this  Work 

The  Hough  transform  has  a prominent  position  in  the  fields  of  image  processing 
and  pattern  recognition  [54,  71].  A great  deal  of  research  has  been  devoted  to  the 
study  of  the  factors  that  contribute  to  the  errors  of  the  Hough  transform  [53,  77,  78, 
88].  A large  number  of  publications  propose  the  use  of  the  Hough  transform  for  the 
task  of  shape  recognition  and  localization  [19,  30].  The  purpose  of  this  work  is  to  add 
to  the  body  of  knowledge  by  studying  the  digitization  errors  of  the  Hough  transform 
with  the  use  of  Digital  Topology  theory.  Digital  Topology  theory  provides  us  with 
the  tools  to  model  digital  straight  line  segments  as  objects  having  a finite  (nonzero) 
width.  The  modeling  of  a digital  straight  line  segment  helps  us  in  the  derivation  of 
a corrected  formula  for  the  worst  case  spreading  of  a digital  straight  line  segment  of 
a given  length.  Other  somewhat  surprising  effects  are  explained  if  one  uses  Digital 
Topology  theory  for  the  study  of  the  digitization  errors.  An  example  is  the  fact  that 
even  in  the  perfect  case  (no  noise,  perfect  digital  straight  line  segments)  the  peak  cell 
may  not  best  represent  the  parameters  of  a digital  straight  line  segment.  The  study 
of  the  performance  prediction  problem  in  the  discrete  image  and  parameter  spaces 
has  been  considered  a very  difficult  problem  to  tackle  [32,  60].  This  is  the  reason  why 
research  efforts  were  directed  towards  the  study  of  this  problem  in  the  continuous 
space  (both  image  and  parameter).  We  show  that  if  one  considers  the  results  of 
Digital  Topology  theory  the  solution  of  this  problem  becomes  straightforward. 

We  show  the  applicability  of  our  results  by  proposing  a new  filtering  technique 
of  the  Hough  space  that  is  based  on  the  estimation  of  the  length  of  the  digital  straight 
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line  segments.  We  also  propose  a new  algorithm  that  uses  the  information  provided 
by  the  Hough  transform  to  recognize  convex  polygons.  Byproducts  of  this  work  are 
the  proof  of  the  uniqueness  of  the  representation  of  a convex  polygon  by  the  peaks  it 
gives  in  the  Hough  space,  and  an  efficient  algorithm  that  constructs  a convex  polygon 
from  the  parameters  of  its  edges. 

It  may  be  thought  that  there  is  too  much  theory  presented  or  that  some  of  the 
theory  is  not  directly  relevant  to  the  applications.  In  our  opinion  an  algorithm  should 
have  a sound  mathematical/theoretical  basis  and  this  is  the  reason  why  we  strove 
for  mathematically/theoretically  sound  algorithms.  A large  number  of  experimental 
results  accompanies  each  algorithm  to  demonstrate  its  soundness.  We  also  strove  for 
a clarity  of  the  ideas  presented  by  explicitly  stating  the  assumptions  underlying  each 
new  theoretical  result.  Care  was  also  taken  so  that  the  parameters  in  the  presented 
experiments  are  explicitly  stated  to  facilitate  the  repeatability  of  them.  On  the 
other  hand  when  we  thought  that  stating  too  many  details  would  not  add  to  the 
understanding,  we  omitted  them. 

1.2  Organization  of  the  material 

This  thesis  is  organized  as  follows.  In  the  remainder  of  this  chapter  we  present  an 
introduction  to  the  Hough  transform,  its  properties  and  review  earlier  works  relevant 
to  this  thesis.  We  also  review  Digital  Topology  theory  which  will  help  us  improve 
our  understanding  of  the  digitization  errors.  The  effects  of  digitization  errors  on 
the  Hough  transform  are  discussed  in  Chapter  2.  In  Chapter  3 we  use  the  results 
of  Chapter  2 to  propose  an  algorithm  that  uses  information  obtained  through  the 
Hough  transform  to  estimate  accurately  the  length  of  a Digital  Straight  Line  (DSL) 
segment.  In  Chapter  4 we  use  the  results  of  the  previous  two  chapters  for  the  recog- 
nition of  convex  polygons.  In  Chapter  5 we  provide  a proof  of  the  uniqueness  of  the 
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representation  of  a convex  polygon  by  the  peaks  it  gives  in  the  Hough  space.  We 
also  propose  an  efficient  algorithm  for  the  reconstruction  of  a convex  polygon  by  the 
equations  of  its  edges.  Discussion  for  the  non-convex  polygon  case  is  also  provided. 
In  Chapter  6 we  present  conclusions  and  future  research  directions. 

1.3  The  Hough  Transform 

1.3.1  Basic  Ideas 

The  Hough  transform  is  a powerful  technique  for  straight  line  detection  in  a 
digital  image.  It  was  first  introduced  as  a patent  by  P.V.C.  Hough  [31].  Rosenfeld 
[74]  was  the  first  to  introduce  it  to  the  image  processing  community. 

The  main  idea  of  the  technique  can  be  understood  if  one  considers  the  problem 
of  the  recognition  of  straight  lines  in  an  image.  The  standard  straight  line  equation: 

y = cx  + m (1.1) 

can  be  rewritten  in  the  format 

m = (-x)c  + y.  (1.2) 

It  can  be  seen  that  each  point  in  the  (x,y)  space  corresponds  to  a straight  line  in 
the  (c,  m)  space.  This  line  represents  the  pencil  of  the  lines  that  pass  through  that 
point.  Obviously,  from  the  same  equation  we  can  conclude  that  each  straight  line  in 
the  ( x,y ) space  is  represented  by  a point  in  the  (c,m)  space.  Practically,  one  cannot 
consider  the  huge  number  of  straight  lines  that  may  appear  in  an  image  [86].  Hence, 
there  is  a need  for  a segmentation  of  the  parameter  space.  An  accumulator  matrix  C is 
initizalized  to  zero.  This  matrix  holds  the  votes  for  the  members  of  the  segmentation 
set  of  the  parameter  space.  Each  feature  point  votes  for  the  lines  (members  of  the 
segmentation  set  of  the  parameter  space)  that  pass  through  it.  Highly  voted  lines 
give  an  indication  of  the  existence  of  a digital  straight  line  segment. 
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Figure  1.1.  Parameters  of  a straight  line  used  in  the  Duda-Hart  parameterization  for 
the  Hough  transform.  Notice  that  for  the  particular  straight  line  9 = </>  - 7r/2. 

Recent  reviews  on  theoretical  developments  and  applications  of  the  Hough  trans- 
form can  be  found  in  the  works  of  Illingworth  and  Kittler  [39]  and  Leavers  [54], 

1.3.2  The  (p.  9)  parametrization 

A disadvantage  of  the  standard  parameterization  of  a straight  line  (Equation 
1.2),  is  that  the  search  space  extends  to  infinity.  This  is  the  reason  why  Duda  and 
Hart  [24]  introduced  the  normal  form  for  the  Hough  transform: 

p = x cos  9 + y sin  9.  (1.3) 

Each  straight  line  is  represented  by  the  distance  p from  the  origin  of  the  Cartesian 
plane,  and  the  angle  9 of  the  perpendicular  to  the  line,  with  the  £-axis  (Figure  1.1). 
The  advantage  of  this  parametrization  is  that  there  is  no  need  for  a search  on  a field 
that  extends  to  infinity.  This  happens  because  9 £ [0, 7r)  and  p £ [-N/y/2,  N/y/2], 
for  an  image  of  N by  N pixels  (see  Figure  1.2). 
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Figure  1.2.  Both  parameters  for  the  (p,9)  parameterization  of  a straight  line  have 
finite  ranges. 

The  Hough  transform  algorithm  that  follows  is  a modified  version  of  the  one 
presented  by  Lam  et  al.  [53]: 
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Algorithm  1 Hough  transform 

select  segmentation/quantization  steps  A p and  A 6 
initialize  C to  zero 
for  all  feature  points  (xm,ym) 
for  j = 0 to  N — 1 
Oj  = i A0 

evaluate  pt  = (xm  — xo)  cos  Oj  + ( ym  — y0 ) sin  03 
Pi  = quantize(pt , A p) 
i = round(pi/ Ap) 

C{f,j)  = C(i,j)  + 1 

end 

end 

end 

where 

quantize(pti  Ap)  — [ ~ ^ ^ ^ ^"-1  < ' 

{ pn  = nAp  otherwise , 

and 

f dn-i  < Pt  < Pn , 

\ Pn  Pn— 1 — Ap. 

With  (a:0,?/o)  we  denote  the  coordinates  of  the  point  of  origin.  A significant  number 
of  alternative  parameterizations  have  been  proposed,  like  the  “Muff  transform”  by 
Wallace  [92]  and  the  “foot  of  the  normal”  by  Davies  [17].  If  one  considers  the  pa- 
rameterizations used  in  the  literature  during  the  last  20  years,  he  would  see  that  the 
overwhelming  majority  of  the  research  community  has  adopted  the  parameterization 
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proposed  by  Duda  and  Hart.  For  an  evaluation  of  the  performance  of  a number  of 
parameterizations  one  should  consult  the  work  of  Hu  and  Ma  [33].  The  properties  of 
the  parameterization  proposed  by  Duda  and  Hart  have  been  summarized  by  Lam  et 
al.  [53]: 


1.  A point  in  the  image  space  corresponds  to  a sinusoid  in  the  parameter  space. 

2.  A point  in  the  parameter  space  corresponds  to  a straight  line  in  the  image 
space. 

3.  Points  lying  on  the  same  straight  line  in  the  image  space  correspond  to  sinusoids 
through  a common  point  in  the  parameter  space. 

4.  Points  lying  on  the  same  sinusoid  in  the  parameter  space  correspond  to  lines 
through  the  same  point  in  the  image  space. 


The  first  property  can  be  proved  as  follows:  for  Ecjuation  1.3  for  the  case  when 
V 7^  0 if  we  substitute  x — y tan/?  we  will  get  p = sin(0  + /?)  which  represents  a 
sinusoidal  curve  in  the  parameter  domain. 

1.3.3  Digitization  Errors 


vanVeen  and  Groen  [88]  were  the  first  to  study  the  spreading  of  the  votes  of  a 
digital  straight  line  segment  to  several  strips.  By  making  the  assumption  that  a DSL 
segment  can  be  modeled  as  a Euclidean  line  segment  (this  means  that  the  centers  of 
the  pixels  that  belong  to  the  DSL  segment  are  on  the  same  straight  line,  Figure  1.3) 
they  derived  the  following  relation: 


n, 


= L 


/ sin(|A0) 

Ap 


J + 2, 


(1.4) 


where  the  notation  [z\  denotes  the  largest  integer  smaller  than  2.  The  symbol  np 
denotes  the  maximum  number  of  strips  a DSL  segment  is  spread  out  (worst  case 
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Figure  1.3.  The  infinitely  thin  Digital  Straight  Line  segment  AB  spreads  its  votes  to 
the  strips  So,  Si,  S2. 

spreading  [64])  and  / is  the  length  of  the  DSL  segment.  Given  the  maximum  length 
of  a line  lmax  one  may  choose  the  Ap  as 

A p > lmax  sin(-A0)  (1.5) 

so  that  the  worst  case  spreading  for  a DSL  segment  of  length  lmax  is  equal  to  2 
(minimum  possible). 

1.3.4  Noise  Modeling 

In  our  effort  to  benchmark  the  algorithms  we  will  rely  on  simulated  images 
where  we  have  control  of  the  parameters  of  the  image.  An  accurate  modeling  of 
the  objects  that  may  exist  in  the  image  and  of  the  various  sources  of  uncertainty 
is  instrumental  for  the  identification  of  the  shortcomings  and  the  advantages  of  the 
proposed  techniques.  This  helps  in  the  understanding  of  the  problem  and  makes 
easier  the  generalization  of  the  techniques  to  deal  with  realistic  situations. 

In  modeling  the  images  that  are  input  to  the  HT  we  had  to  take  into  account 
two  sources  of  uncertainty,  namely: 

1.  noise , noise  pixels  and  pixels  belonging  to  objects  other  than  the  object  of 
interest  (in  our  case,  other  than  the  digital  straight  line  under  consideration), 
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2.  straight  line  digitization , the  straight  lines  detected  using  conventional  edge 
detectors  are  not  perfect  and  do  not  satisfy  the  conditions  imposed  by  Digital 
Topology  theory  (see  Section  1.5).  These  line  segments  often  appear  as  broken 
segments,  with  changing  directions  and  often  they  are  rounded  at  the  ends 
deviating  from  the  straight  line  model  [18]. 

The  model  we  chose  for  the  noise  in  the  binary,  edge  detected  images  was  that 
of  a uniform  distribution.  This  is  a very  reasonable  choice  since  objects  appear  in 
all  parts  of  the  image  with  equal  probability.  We  also  point  out  that  in  most  of 
the  published  work  that  used  synthetic  images  to  benchmark  the  performance  of 
various  implementations  of  the  Hough  transform  the  noise  was  uniformly  distributed 
[4,  9,  13,  16,  34,  38,  52,  56,  60].  In  this  work  when  we  refer  to  a p uniform  noise  we 
mean  that  the  probability  of  a pixel  to  be  noisy  is  equal  to  100p%.  More  discussion 
on  this  will  appear  in  Section  2.6. 

Another  noise  model  that  has  found  some  use  is  the  one  studied  by  Hu  and  Ma 
[34].  This  noise  model  assumed  a Gaussian  distribution  of  the  background  noise. 

1.3.5  Postprocessing  of  the  Hough  Transform 

The  evaluation  of  the  information  provided  by  the  Hough  transform  is  not  an 
easy  task.  Several  methods  have  been  proposed  for  the  accurate  detection  of  maxima 
in  the  parameter  space.  Leavers  and  Boyce  [55]  have  put  forward  the  idea  of  filtering 
of  the  Hough  space  before  detecting  the  maxima.  Their  filter  derivation  incorporated 
knowledge  about  the  dimensions  of  the  features  under  detection  (i.e.,  the  length  of 
the  digital  straight  line  segment)  and  the  quantization  of  both  the  image  and  the 
parameter  space.  They  showed  that  the  maxima  in  Hough  space  corresponding  to 
straight  line  segments  give  butterfly  features  (see  Figure  1.4).  Using  this  information 
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they  derived  the  following  filter: 


t = 


0-2  0 
1 2 1 
0-2  0 


(1.6) 


for  improved  detectability  and  suppression  of  noise.  Princen  et  al.  [71]  considered  the 
design  of  optimum  ID  filters  towards  the  above  mentioned  goal.  In  another  paper 
[70]  the  same  authors  proposed  a sliding  along  the  p-direction  of  the  filter  (as  in 
Figure  1.5): 


t = 


(1.7) 


to  suppress  local  peaks  due  to  noisy  features. 

An  interesting  idea  has  been  proposed  by  van  Veen  and  Groen  [88].  By  assuming 
a known  length  they  used  Equation  1.4  to  calculate  np.  They  proposed  the  use  of 
a filter  of  length  np  that  takes  into  account  the  votes  that  the  digital  straight  line 
segments  spreads  to  a set  of  parallel  strips  to  be  slid  along  the  p direction.  For 
example  when  np  — 3;  the  proposed  filter  is  the  following: 


t = 


1 

1 

1 


(1.8) 


A big  problem  of  these  filtering  techniques  is  that,  after  filtering,  thresholding  is 
needed  and  there  is  no  intuitive  approach  towards  the  selection  of  the  threshold  value. 
Another  problem  is  that  these  filters  incorporate  knowledge  about  the  dimensions  of 
the  features  under  detection.  This  kind  of  information  is  not  always  available. 


1.4  Digital  Topology 

Digital  Topology  provides  the  foundations  for  the  study  of  the  topological  prop- 
erties of  digital  images.  The  results  from  this  theory  will  be  particularly  useful  in  our 
study  of  the  effects  of  the  digitization  errors  on  the  Hough  transform.  The  branch 
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(a) 


(b) 

Figure  1.4.  Image  containing  a single  DSL  segment  and  the  butterfly  formation 
it  gives  in  the  Hough  domain.  The  size  of  the  image  is  256  by  256  pixels.  The 
discretization  steps  of  choice  were  A p = h and  A 9 = 7t/256,  where  h is  the  size  of 
the  pixel. 
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Figure  1.5.  Illustration  of  the  filtering  with  the  template  proposed  by  Princen  et  al. 
The  distances  of  the  parallel  strips  »Sb,  S\ , and  S2  from  the  point  of  origin  O are 
Pi  + A p and  pi  + 2A p. 

of  Digital  Topology  that  deals  with  unsegmented  gray  scale  images  (Fuzzy  Digital 
Topology  [51])  will  not  be  discussed  in  this  section.  We  will  only  concentrate  on  re- 
sults that  are  relevant  to  the  definition  and  modeling  of  digital  straight  line  segments 
in  binary  images. 

The  need  for  the  use  of  Digital  Topology  theory  for  the  study  of  the  digitiza- 
tion errors  is  a natural  one,  if  one  considers  that  the  Hough  transform  is  actually 
a template  matching  technique  [82,  85].  Therefore,  there  is  a need  for  an  accurate 
modeling  of  the  object  (straight  lines  in  this  case)  one  is  trying  to  match  (detect). 
This  modeling  is  provided  by  the  results  of  Digital  Topology  theory. 

The  digitization  scheme  used  in  this  work  is  the  standard  grid  intersection 
method  [28]:  given  a coordinate  grid  superimposed  on  the  curve,  then  whenever 
the  curve  crosses  a grid  line,  the  grid  point  nearest  to  the  crossing  becomes  a point 
of  the  digitization  of  the  curve.  We  will  assume  that  the  size  of  a pixel  is  h by  h.  It 
has  been  shown  [21]  that  for  straight  lines  the  two  quantization  methods,  the  grid 
intersection  and  the  object  boundary  quantization,  are  equivalent.  According  to  the 
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Figure  1.6.  The  8 basic  directions  of  the  image.  Points  (0,2,4, 6)  form  the  4- 
neighborhood  of  0 while  points  (0  — 7)  form  the  8-neighborhood  of  that  point. 

object  boundary  quantization,  the  outermost  points  of  an  object  are  digitized  by  a 
chaincode  string  [21]. 

We,  now,  introduce  the  important  idea  of  adjacency  in  the  two-dimensional 
image  space. 

Definition  1.1  [51]  Two  lattice  points  in  the  plane  are  said  to  be  8-adjacent  if  they 
are  distinct  and  each  coordinate  of  one  differs  from  the  corresponding  coordinate  of 
the  other  by  at  most  1. 

Definition  1.2  [51]  Two  lattice  points  in  the  plane  are  said  to  be  A-adjacent  if  they 
are  8—  adjacent  and  differ  in  at  most  one  of  their  coordinates. 

Figure  1.6  illustrates  the  8 and  4-adjacency.  In  this  figure,  points  in  the  direc- 
tions 0,2,4  and  6 are  4-adjacent  to  point  O.  These  points  along  with  the  points  in 
the  directions  1,3,5  and  7 constitute  the  8-neighborhood  of  point  O.  Directions  0-7 
are  the  8 basic  directions  of  the  digital  image. 

Throughout  this  work  we  will  be  dealing  with  binary  images.  These  images 
are  the  result  of  the  edge  detection-thresholding  of  gray  scale  images.  The  whole 
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process  of  edge  detection  and  thresholding  is  used  for  the  demarkation  of  the  pixels 
to  feature  pixels  (i.e. , pixels  where  the  change  of  intensity  is  large)  and  non-feature 
pixels  (pixels  where  no  big  changes  take  place).  Digital  Topology  theorists  took  a 
different  approach  when  defining  the  feature  points  than  the  members  of  the  image 
processing  community.  The  former  regarded  the  black  points  as  the  feature  points 
while  the  latter  used  the  white  points  as  feature  points.  In  the  course  of  this  work, 
we  will  be  complying  with  the  latter. 

Definition  1.3  [51]  A lattice  point  that  has  a value  of  1 is  called  a white  (feature) 
point. 

The  difference  between  this  definition  and  the  one  provided  by  Kong  and  Rosenfeld 
[51]  is  that  according  to  our  scheme,  feature  points  are  white  points  while  for  the 
reference  mentioned  above,  black  points  are  considered  feature  points.  We  now  define 
the  property  of  connectivity  for  sets  of  feature  points.  This  property  is  fundamental 
for  the  definition  of  a digital  image. 

Definition  l.f  [51]  It  can  be  said  a set  S of  feature  points  is  8-connected  if  S cannot 
be  partitioned  into  two  subsets  that  are  not  8-adjacent  to  each  other. 

An  example  of  an  8-connected  set  of  feature  points  is  shown  in  Figure  1.7.  We 
now  provide  a definition  of  a digital  image  as  it  will  be  used  in  this  work.  For  a more 
general  definition  one  can  refer  to  [51]. 

Definition  1.5  [51]  A digital  image  is  a quadruple  (V,8,4,B)  where  V = Z2 , B C Z2 . 

Z is  the  set  of  integers.  This  definition  says  that  the  8— adjacency  rule  will  be  used 
for  white  components  while  the  4— adjacency  will  be  used  for  black  components.  The 
need  for  a dual  adjacency  rule  was  put  forward  to  tackle  a paradox  first  discussed  by 
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Figure  1.7.  An  8-connected  set.  Point  C is  a non-simple  point  while  points  A and  B 
are  simple  points.  If  C is  deleted  the  numbers  of  feature  and  non-feature  components 
remain  the  same.  If  point  B is  deleted  we  get  two  feature  components.  The  first  is 
point  C and  the  second  the  rest  of  the  object.  The  same  happens  if  we  delete  point 
A. 
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Duda  et  al.  [51].  The  definitions  that  follow  will  be  used  in  the  definition  of  a digital 
straight  line  segment. 

Definition  1.6  [51]  A point  P is  a simple  point  of  a digital  picture  if  and  only  if  both 
the  number  of  feature  components  and  the  number  of  non-feature  components  stay  the 
same  when  P is  deleted. 

Definition  1.7  [51]  A simple  digital  arc  is  a connected  set  of  white  (feature)  points 
each  of  which  is  adjacent  to  just  two  feature  points  with  the  exception  of  two,  the  end 
points  of  the  arc,  that  are  each  adjacent  to  just  one  other  feature  point. 

In  Figure  1.7  point  C is  a simple  point  while  points  A,  B are  not  simple  points. 
Figure  1.8  shows  an  example  of  a simple  digital  arc. 

A popular  representation  of  a simple  digital  arc  is  the  chain  code.  The  chain  code 
specifies  the  direction  of  the  contour  at  each  pixel  of  the  digital  arc  [43].  Directions 
are  quantized  into  the  set  of  the  eight  basic  directions  (see  Figure  1.6).  For  example, 
the  chain  code  of  the  digital  arc  shown  in  Figure  1.8,  starting  from  the  endpoint  Cfi 
is  6766770001.  The  chain  code  specifies  a digital  arc  up  to  a translation  factor.  If 
one  wants  to  determine  the  position  of  the  digital  arc  in  the  image  space  he  also  has 
to  know  the  coordinates  of  one  of  its  endpoints. 

Definition  1.8  [j8]  A digital  arc  R is  said  to  be  a digital  straight  line  if  there  exists 
a directed  line  segment  f whose  digital  image  is  R. 

The  following  lemma  will  be  used  throughout  this  thesis  for  the  modeling  of 
digital  straight  line  segments  and  digital  straight  lines.  This  modeling  will  be  very 
useful  for  the  study  of  the  digitization  errors  of  the  Hough  transform. 

Lemma  1.1  . [f9]  Let  (s,6)  and  (s  + As,6)  be  two  parallel  lines  such  that  the  smaller 
of  the  two  distances,  the  horizontal  and  the  vertical,  between  them  is  equal  to  h.  S 
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Figure  1.8.  A digital  arc.  Notice  that  except  for  the  two  endpoints  0\  and  02  all 
other  points  have  two  points  of  their  8-neighborhood  belonging  to  the  digital  arc. 
Points  0\  and  02  have  only  one  point  of  their  8-neighborhood  belonging  to  the  digital 
arc. 
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Figure  1.9.  According  to  Digital  Topology  theory  a DSL  has  a certain  width.  It  is 
obvious  that  the  centers  of  the  pixels  of  the  DSL  segment  between  u and  v do  not 
belong  to  the  same  Euclidean  straight  line. 

is  the  set  of  lattice  points  that  lie  between  the  two  lines  and  on  (s,0)  but  not  on 
(s  + As,9).  If  u and  w are  two  points  of  S and  L is  the  set  of  points  that  belong  to  S 
and  lie  between  u and  w,  u and  w included,  then  L is  a digital  straight  line  segment. 

Figure  1.9  illustrates  the  ideas  of  Lemma  1.1. 

Since  throughout  this  work  we  will  be  dealing  with  length  estimation  of  digital 
straight  line  segments,  a final  word  has  to  be  said  about  the  definition  of  it.  In 
this  work,  “accurate”  length  of  a DSL  segment  will  refer  to  the  distance  between 
its  two  endpoints.  This  estimator  in  literature  [22]  has  been  referred  to  as  the  Most 
Probable  Original  (MPO)  estimator.  Dorst  and  Smeulders  [22]  showed  that  even 
though  the  MPO  estimator  is  simpler  to  compute  than  the  optimum  Best  Linear 
Unbiased  Estimator  it  is  almost  as  accurate. 


CHAPTER  2 

DIGITIZATION  ERRORS  OF  THE  HOUGH  TRANSFORM 

2.1  Introduction 

In  this  chapter  we  will  be  dealing  with  the  subject  of  digitization  errors  of  the 
Hough  transform.  In  particular,  we  provide  an  indepth  study  of  the  errors  due  to  the 
digitization  of  the  image  and  the  discretization  of  the  parameter  space.  Early  studies 
did  not  take  into  account  the  discretization  of  the  image  space  when  they  developed 
their  theory  of  line  segment  spreading  to  several  cells  of  the  parameter  space  [88]. 

Our  treatment  of  digitization  errors  relies  on  geometric  arguments.  Members  of 
the  parameter  space  are  modeled  as  strips  in  the  image  space  and  digital  straight  lines 
are  modeled  as  strips  according  to  Lemma  1.1.  In  Section  2.2  we  treat  the  problem 
of  the  spreading  of  a digital  straight  line  segment.  In  Section  2.3,  an  equation  for 
the  upper  bound  of  the  votes  a strip  takes  due  to  a DSL  segment  is  derived.  In 
Section  2.4  we  use  the  results  of  Section  2.3  to  treat  the  problem  of  peak  extension 
[53].  We  may  say  that  Sections  2.2,  2.3  and  2.4  provide  a theoretical  framework  for 
the  modeling  of  peak  formation.  Section  2.5  contains  a discussion  on  the  choice  of 
A/9.  Many  parameterizations  have  been  proposed,  all  of  them  targeted  at  designing 
an  optimum  A p according  to  criteria  such  as  uniformity  of  noise  in  the  parameter 
space,  better  accuracy,  etc.  [13,  15,  64].  In  the  same  section  we  show  that  according 
to  Digital  Topological  criteria  there  should  be  an  adaptive  A p that  is  a function  of 
the  orientation  of  the  strips.  In  Section  2.6  we  deal  with  the  effects  of  uniform  noise 
on  the  parameter  space. 
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2.2  The  Spreading  of  a DSL  Segment 

2.2.1  Theoretical  Observations 

In  our  effort  to  treat  the  problem  of  the  spreading  of  a DSL  segment,  we  will  use 
results  from  the  Digital  Topology  theory  (see  Chapter  1).  First  we  give  a definition 
of  the  spreading  of  a DSL  segment  to  several  cells. 

Definition  2.1  Consider  a DSL  segment  with  an  angle  6 in  an  image  and  a (Ap,  A0) 
segmentation  of  the  Hough  space.  For  the  j-th  column  of  the  quantization  set,  spread- 
ing is  the  phenomenon  of  the  DSL  segment  splitting  its  votes  to  several  cells  of  this 
column. 

This  definition  of  the  spreading  phenomenon  applies  to  any  column  of  the 
paremeter  set.  With  we  will  denote  the  maximum  possible  cells  a DSL  seg- 
ment spreads  out  for  the  j-th  member  of  the  0-sampling  set  (j-th  column  of  the 
accumulation  matrix). 

Since  we  are  dealing  with  a DSL  segment,  longitudinal  information  is  not  enough 
to  characterize  it.  If  the  object  we  are  dealing  with  was  an  infinitely  thin  straight  line 
segment,  its  length  would  have  been  sufficient  for  the  calculation  of  n^.  From  Lemma 

1.1  we  know  that  a DSL  segment  is  an  object  that  has  a nonzero  width.  If  h is  the  size 
of  a pixel  we  can  say  that  for  each  DSL  segment  with  an  angle  6 we  can  always  find 
a pair  of  parallel  straight  lines  whose  distance  is  equal  to  h max( | cos  0|,  |sin0|)  and 
contains  all  the  pixels  that  belong  to  the  segment  (see  Lemma  1.1).  The  uniqueness 
of  this  pair  of  straight  lines  of  angle  9 is  not  guaranteed  though.  It  can  be  proven 
that  for  a DSL  segment  of  angle  9 there  can  be  found  an  infinite  number  of  pairs  of 
straight  lines  that  define  a DSL  and  contain  the  segment  (see  Lemma  1.1).  Another 
definition  is  given: 
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Definition  2.2  The  member  9j  = jA9  of  the  angle  quantization  set  which  is  closest  to 
the  angle  9S  of  the  DSL  segment  will  be  called  the  primary  angle  of  the  DSL  segment 
for  the  given  segmentation.  The  strips  of  angle  90  that  get  voted  by  the  pixels  of  the 
DSL  segment  form  the  set  of  primary  cells/strips  of  the  DSL  segment  for  the  given 
segmentation. 


It  is  to  be  noticed  that  we  did  not  define  the  set  of  the  primary  strips  either  as 
the  set  of  the  strips  whose  member  holds  the  maximum  number  of  votes  due  to  the 
segment  or  as  the  set  of  strips  where  the  minimum  n ^ appears.  We  took  this  approach 
because,  as  we  will  prove  later,  it  is  not  guaranteed  that  the  peak  cell  (highest  voted 
cell)  belongs  to  the  set  of  primary  strips  or  that  the  set  of  primary  strips  has  the 
fewest  number  of  members  (smallest  nJp). 

We  now  proceed  to  derive  a formula  for  the  worst  case  spreading  of  a DSL 
segment,  np.  First  we  define  the  spreading  width. 


Definition  2.3  Spreading  width,  w,  of  a DSL  segment  is  the  length  of  the  projection 
of  the  segment  to  an  axis  that  is  orthogonal  to  the  set  of  primary  strips. 


The  sources  that  contribute  to  the  spreading  of  a DSL  segment  are  two:  the 
length  and  the  width  of  the  segment.  According  to  what  we  said,  the  width  of  the 
segment  can  be  set  equal  to  h max( | cos  9S | , |sin#s|).  If  we  were  talking  about  a Eu- 
clidean straight  line  segment,  the  maximum  possible  spreading  width  would  have  been 
/sin(|A0)  where  / is  the  length  of  the  segment  (see  Figure  1.3).  For  digital  straight 
lines  the  spreading  width  of  the  digital  segment  is  bounded  by  the  sum  of  the  projec- 
tion of  its  length  /sin(|A0)  and  the  projection  of  the  width  hmax(|  cos^s|,  | sin0s|)- 
cos(^A$)  (Figure  2.1).  Therefore,  the  maximum  possible  value  the  spreading  width 
can  take  is  given  by 


w-max  — lsm(-A9)  + hmax(|cos0s|,  | sin  0S|)  cos(-A0). 


(2.1) 


22 


hmax( I cos9s I , I sin0g I ) / 2 


Figure  2.1.  The  spreading  width,  w,  of  a DSL  segment  to  a set  of  strips  belonging 
to  the  same  column  of  the  parameter  space  is  bounded  by  the  sum  of  the  projection 
of  the  length  of  the  segment  lsm(86)  and  the  projection  of  the  width  of  the  seg- 
ment h max(|  cos  0S|,  | sin  0S|)  cos(<50).  For  the  DSL  segment  of  this  figure,  / = \AB\, 
l max( | cos  6S | , | sin 6S | ) cos  60  = \APBP\  and  w = \AxBi\. 


This  leads  to  the  following  equation: 


w „ 


n, 


_ | _max  | 2 _ I 

L Ap  J + L A p 


l sin(|A0)  /imax(|cos  0a|,  | sin  ^5|)  cos(|A6*) 


+ 


A P 


J+2.  (2.2) 


If  we  recognize  that  max(cos(§A0)  max(|  cos0,|,  | sin0s|))  = 1 we  will  get 

/sin(iA0)  h 


(2.3) 


If  we  compare  Equation  2.3  with  Equation  1.4  of  van  Veen  and  Groen  we  can  see  that 
there  is  one  additional  term  within  the  brackets  which  accounts  for  the  width  of  the 
digital  straight  line  segment. 

2.2.2  Experimental  Results 


Several  experiments  were  prerformed  to  prove  the  claim  that  Equation  2.3  rep- 
resents an  improvement  over  Equation  1.4.  Consider  the  DSL  segment  with  vertices: 


23 


(43, 13)  — (204, 24).  The  chaincode  of  this  DSL  segment  is  00000001000000000000010- 
0000000000000100000000000000100000000000001000000000000001000000000000001- 
00000000000001000000000000001000000000000001000000000000010000000.  A com- 
parison of  the  prediction  of  Equation  2.3  and  Equation  1.4  with  varying  A p is  pre- 
sented in  Table  2.1.  Obviously  as  A p becomes  smaller  the  second  term  of  Equation 
2.3  becomes  more  important  and  the  discrepancy  between  Equation  1.4  and  the  mea- 
sured values  becomes  larger. 

Table  2.1.  Predictions  of  the  worst  case  spreading  of  the  model  of  van  Veen  and 
Groen  with  the  model  presented  in  this  section  for  a number  of  Ap  segmentations. 
The  fourth  column  reports  the  measured  spreading.  A 6 is  equal  to  7t/256. 


Ap 

van  Veen  and  Groen 

New  Model 

Measured  Spreading 

h 

2 

3 

3 

h/2 

3 

5 

4 

h/ 4 

5 

9 

8 

h/6 

7 

13 

11 

h/8 

9 

17 

14 

2.3  A Model  for  the  Prediction  of  the  Votes  of  a Strip  Due  to  a DSL  Segment 
2.3.1  Theory 

The  analysis  that  follows  will  give  a theoretical  prediction  for  the  number  of 
votes  a strip  takes  due  to  a DSL  segment  as  a function  of  the  angle  between  them. 
Suppose  that  the  angle  of  the  segment  is  #s,  and  the  angle  of  the  strip  is  0j  = jA6. 
The  angle  between  the  segment  and  the  strip  is  86  = 1 9j  — 6S\.  From  Figure  2.2  we 
can  see  that  the  area  of  the  intersection  of  a line  L and  a strip  S is  equal  to 

^ (Ap)(h)max(\ cos0s|,  | sint^l) 
tan(^0) 

From  Lemma  1.1,  for  any  strip  with  an  arbitrary  angle  6 for  a width  equal  to 
h max(|  cos  6\,  | sin  #|)  corresponds  1 pixel  per  h/  max(|  cos  0|,  | sin  6\)  length  (by  using 
the  term  pixel  we  mean  center  of  the  pixel).  This  is  the  only  case  where  pixels  are 
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hmax ( I cos (0) I , I sin (0) | ) 


Figure  2.2.  Area  of  intersection  of  a line  (L)  and  a strip  (S).  Their  relative  angle  is 
equal  to  86. 


distributed  uniformly  among  parallel  strips  (uniform  distribution  of  the  pixels  means 
constant  number  of  pixels  per  unit  area).  Even  though  this  is  the  case  we  will  assume 
that  the  distribution  of  pixels  among  parallel  strips  is  fairly  uniform,  in  other  words 
there  are  no  big  differences  in  the  number  of  points  between  parallel  strips  of  the 
same  length.  This  means  that  for  all  strips  there  is  (1  pixel) /h2  uniformly. 

If  we  use  of  the  assumption  discussed  previously,  we  find  that  the  maximum 
possible  number  of  pixels  of  a DSL  that  belong  to  a strip  with  A p — h is  equal  to: 


Vmax 


max(|  cos  6S\,  | sin  0S|) 
tan(<5$) 


(2.5) 


A simple  example  showing  that  our  assumption  about  the  uniformity  of  the 
distribution  of  pixels  does  not  always  hold  is  shown  in  Figure  2.3.  In  this  figure  it  is 
shown  that  two  objects  that  have  the  same  area  contain  different  number  of  pixels. 


2.3.2  Experimental  Results 


In  Figure  2.4,  plots  are  shown  of  the  maximum  number  of  votes  of  a strip 
vmax  as  a function  of  its  angle  relative  to  the  primary  strip  for  three  DSL  segments 
with  angles:  (07T,  207t/256,  647t/256).  The  Hough  transform  of  an  image  containing  a 
single  digital  straight  line  segment  was  performed  and  the  highest  voted  cells  for  each 
column  (relative  angle  with  the  segment)  were  plotted.  The  A p/h  ratio  was  equal  to 
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Figure  2.3.  Illustration  of  the  fact  that  if  a triangle  is  translated  it  will  not  contain 
the  same  number  of  pixels.  Triangles  A and  B have  the  same  area  but  contain  a 
different  number  of  pixels. 

one.  According  to  the  predictions  of  Equation  2.5  the  number  of  votes  has  an  inverse 
tangential  dependence  on  the  relative  angle  60  of  the  line  with  the  strip. 

The  study  of  the  distribution  of  the  votes  of  a segment  into  parallel  strips  (spread 
out  of  the  segment)  proved  to  be  very  useful  for  two  reasons.  First,  it  provides  an 
insight  of  the  spreading  of  a DSL  segment  as  it  gets  longer.  Second  it  helps  us  confirm 
experimentally  the  theoretical  observations  of  this  subsection.  Figure  2.5(a)  shows 
how  the  votes  of  a segment  get  distributed  into  four  different  strips  as  it  becomes 
longer.  In  the  course  of  study  of  strip  C of  this  figure  it  can  be  observed  that  the 
plot  can  be  divided  into  five  different  regions.  In  region  I (0-50)  strip  C does  not  get 
voted.  In  region  11(50-130)  it  shares  the  votes  coming  from  the  segment  with  strip  B. 
Strips  B and  C have  the  same  orientation  and  their  distances  from  the  origin  differ 
by  h.  In  region  III  (130-140)  the  C strip  is  the  exclusive  recipient  of  the  votes,  while 
in  region  IV  (140-180)  strips  C and  D share  the  votes.  Finally  in  region  V no  more 
votes  are  added  to  strip  C.  The  most  important  thing  to  notice  is  that  in  regions  II 
and  IV  the  new  votes  go  alternatively  to  B and  C (for  region  II),  and  C and  D (for 
region  IV). 
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(a) 


(b) 


(c) 

Figure  2.4.  Votes  of  strips  as  a function  of  their  angle  relatively  to  the  DSL  for: 
(a)  a segment  of  07r  orientation  (b)  a segment  of  207r/256  orientation  (c)  a segment 
of  647t/256  orientation.  The  continuous  line  indicates  the  theoretical  predictions 
(Equation  2.5). 
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Special  cases  appear  for  Oj  = 0 and  6j  - 7t/4  (see  Figures  2.5(b)-(c)).  In  these 
cases  there  are  only  three  regions.  These  special  cases  arise  because  the  votes  cannot 
alternate  between  two  neighboring  strips  and  it  is  the  result  of  the  two  directions 
being  coincident  with  the  two  of  the  basic  directions  of  the  image  (Figure  1.6). 

2.4  Peak  Extension 

2.4.1  Theory 

Peak  extension  is  the  effect  of  the  propagation  of  the  peak  to  the  ^-direction  of 
the  parameter  space.  Two  strips  having  small  relative  angle  share  a large  number  of 
pixels.  Therefore,  if  one  of  them  gets  a high  number  of  votes  due  to  a straight  line 
segment  the  other  may  also  get  highly  voted.  The  purpose  of  this  subsection  is  to 
study  the  peak  extension  and  provide  a new  quantitative  and  qualitative  insight. 

A definition  similar  to  the  one  given  in  Subsection  2.2  is  given: 

Definition  2-4  The  set  of  strips  whose  relative  angle  with  the  members  of  the  set  of 
primary  strips  is  equal  to  ± A 9 will  be  called  the  set  of  secondary  strips. 

Hence,  the  members  of  the  accumulation  matrix  that  holds  the  votes  of  the 
secondary  strips  belong  to  columns  whose  index  differs  from  the  index  of  the  column 
of  the  set  of  primary  strips  by  one  (their  angle  is  equal  to  9j  ± AO). 

Even  though  the  spreading  of  a DSL  segment  is  almost  unavoidable,  let  us 
assume  for  the  moment  that  a DSL  segment  belongs  to  a single  strip  (p,-,0j).  It  can 
easily  be  seen  from  Figure  2.6  that  cases  may  exist  that  even  though  a DSL  segment 
belongs  totally  to  a strip  (i.e.,  it  is  exactly  collinear  and  has  exactly  the  same  p with 
it),  the  number  of  votes  given  to  that  strip  is  equal  to  the  number  of  votes  given  to 
a secondary  strip.  The  maximum  number  of  votes  a strip  can  take  due  to  a DSL 
segment  is  given  by  Equation  2.5.  The  fact  that  the  segment  is  exactly  collinear  with 
the  strip  ( pi,0j ) means  that  the  relative  angle  of  the  DSL  segment  with  the  strips  is 
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(a) 


(b) 


(c) 

Figure  2.5.  Number  of  votes  primary  strips  get  as  a function  of  the  number  of 
pixels  belonging  to  a digital  straight  line  segment:  (a)  for  spreading  out  into  strips 
of  orientation  207r/256.  Five  regions  can  be  distinguished;  (b)  for  spreading  out 
into  strips  having  orientation  07t/256.  Only  three  regions  can  be  distinguished;  (c) 
for  spreading  out  into  strips  having  orientation  tt/4.  Only  three  regions  can  be 
distinguished.  This  Figure  also  illustrates  the  existence  of  the  strips  with  orientation 
7t/4  having  simple  points  throughout  their  length.  Strip  B gets  double  the  number 
of  votes  predicted  by  the  theory  presented  in  Section  2.3.  This  happens  because 
the  digital  straight  line  segment,  whose  votes  it  shares,  is  partitioned  among  digital 
straight  line  segments  having  the  basic  direction  7r/4.  Two  of  those  happen  to  belong 
to  the  strip  B and  this  is  the  reason  why  it  takes  double  the  number  of  votes  predicted 
by  the  theory. 
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equal  to  A 8.  Therefore,  a theoretical  limit  for  the  number  of  votes  of  the  secondary 
strips  can  be  found  if  we  use  Equation  2.5: 

_ max(  | cos  | , | sin  8S  | ) 

max  tan(A0)  ‘ (2'6) 

If  a segment  of  angle  8S  has  a number  of  votes  less  than  this  value,  there  may  be 
cases  where  more  than  one  best  estimate  of  (p,  8)  is  given  by  the  HT. 

Another  question  is  under  what  conditions  the  maximum  may  be  attributed  to 
a secondary  strip  instead  of  a primary  strip.  This  is  a very  important  issue  to  address 
because  many  HT  postprocessing  algorithms  are  based  on  the  assumption  that  under 
ideal  circumstances  the  highest  voted  cell  gives  the  best  estimate  of  the  parameters 
of  the  DSL  segment.  In  other  words  the  question  we  need  to  answer  is  the  following: 
is  there  any  spatial  configuration  under  which  the  interplay  of  the  spreading  of  the 
peak  and  the  peak  extension  may  lead  to  a secondary  strip  getting  more  votes  than 
the  highest  voted  member  of  the  primary  set? 

Suppose  that  a DSL  segment  has  an  angle  equal  to  one  of  the  angles  of  the  8- 
sampling  set.  In  this  case  spreading  out  exists  but  if  we  assume  that  a step  A p = h 
has  been  chosen,  the  DSL  segment  can  split  its  votes  between  only  two  neighboring 
strips  (see  Figure  2.?).  The  maximum  number  of  votes  a secondary  strip  can  get  for 
the  case  of  the  DSL  segment  being  collinear  with  the  primary  strips  can  be  calculated 
from  Equation  2.6.  If  a segment  is  spread  into  two  different  primary  strips  and  if 
C\  and  C2  are  the  votes  of  these  strips  for  the  case  where  C\  = C2  — vmax  — 1,  a 
segment  with  as  many  votes  as  2vmax  — 2 can  give  a peak  in  the  secondary  set  of 
strips.  The  conclusion  is  that  even  in  the  ideal  case  (i.e.  no  noise,  perfect 
DSL  segment),  it  is  not  guaranteed  that  the  angle  of  the  peak  is  the  best 
estimate  of  the  angle  of  the  DSL  segment.  In  other  words  the  inequality 


(2.7) 
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Figure  2.6.  Even  though  the  DSL  segment  L is  exactly  collinear  with  strip  Si  the 
votes  of  strips  Si  and  S2  due  to  that  segment  are  equal. 

where  6j  is  the  angle  where  the  peak  occurs  and  0S  is  the  angle  of  the  DSL  segment, 
does  not  necessarily  hold. 

2.4.2  Experimental  Results 

The  DSL  segment  (50, 80)-(  160, 108)  whose  chain  code  is:  010001000100010001- 
0001000100010001000100010001000100010010001000100010001000100010001000100- 
0100010001000100010,  gave  a peak  (highest  voted  cell)  at  (102,150).  The  angle  of 
orientation  for  this  segment,  as  calculated  by  the  coordinates  of  its  endpoints,  is 
0S  = 1.8200.  The  error  in  the  estimation  of  the  orientation  is:  e(0)  > 0.084  > 
0.061  = A0/2. 

2.5  On  the  Choice  of  Ap 

The  effects  of  the  choice  of  the  segmentation  steps  Ap  and  A0  have  been  thor- 
oughly studied  during  the  past  twenty  years  [39].  In  this  section  we  will  treat  the 
choice  of  Ap  from  a Digital  Topological  perspective. 

A sufficient  condition  for  a strip  to  represent  a DSL  is  that  Ap  = /imax(|  cos0|,- 
| sin#|).  This  condition,  provided  by  Lemma  1.1,  is  sufficient  but  not  necessary  for  a 
strip  to  represent  a DSL.  Therefore  an  object  which  is  enclosed  between  two  parallel 
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Figure  2.7.  Two  cases  for  a line  whose  orientation  exactly  coincides  with  the  orien- 
tation of  a member  of  the  parameterization  set.  (a)  Line  L belongs  to  strip  S2-  (b) 
Line  L spreads  its  votes  to  strips  S\  and  S2. 

straight  lines  whose  distance  is  greater /smaller  than  hmax(|  cos#|,  | sin#|)  may  be 
a DSL  segment  but  this  is  not  guaranteed.  If  the  object  enclosed  between  these 
parallel  straight  lines  is  a DSL,  then  from  Lemma  1.1,  it  is  guaranteed  that  we  can 
find  a pair  of  straight  lines  whose  distance  is  h max(|  cos  0|,  | sin  0|)  and  which  contain 
the  DSL  segment  and  define  a DSL.  If  Ap  is  chosen  such  that  A p > hmax(|  cos0|,- 
| sin  ^| ) non-simple  points  may  appear  in  the  digital  object  defined  by  the  strip.  The 
appearance  of  non-simple  points  in  a strip  means  that  in  the  (x,y)  space,  for  an 
angle  9 £ [0,  7t/4),  there  are  columns  where  two  pixels  belong  to  this  strip  as  in 
Figure  2.8.  For  9 £ [7r/4,  7t/2)  the  same  phenomena  appear  but  instead  of  columns 
we  have  rows.  Two  cases  are  of  special  interest  and  will  help  us  obtain  a deeper 
insight.  Suppose  that  we  are  dealing  with  a strip  parallel  to  the  x— axis.  Suppose 
also  that  Ap  = 1.5 h > h max(|  cos  0|,  | sin  #|)  = h (since  9j  = 0).  From  what  has 
already  been  said  it  is  obvious  that  a strip  may  exist  which  has  a simple  point.  For 
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Figure  2.8.  Illustration  of  the  fact  that  if  a strip  L has  a width  larger  than 
hmax(|  cos0|,  | sin0|)  it  may  have  non-simple  points.  Both  pixels  A and  B (their 
centers)  belong  to  the  digital  object  that  is  bounded  by  the  two  parallel  straight 
lines. 

this  case,  the  simple  points  appear  throughout  the  length  of  the  strip  (see  Figure 
2.9). 

If  we  consider  the  case  where  9 = 7t/4  (which  is  one  of  the  basic  directions  too), 
for  a choice  of  A p = 1.5 h > 2hmax(|  cos  0|,  | sin  0|)  strips  that  represent  three  DSLs 
of  angle  9 = 7r/4  and  two  DSLs  of  the  same  slope  would  alternate  (see  Figure  2.10). 
For  the  case  where  A p belongs  to  the  interval  (hmax(|  cos  0\,  | sin  <9|),  2h  max(|  cos0|,- 
| sin  ^|))  phenomena  similar  to  the  ones  discussed  for  the  9 — 0 and  A p — 1.5 h will 
appear.  For  9 6 (0,tt/4)  and  A p £ (h  max(|  cos  6»|,  | sin  0|),  2h  max(|  cos  0\,  | sin  6>|)) 
simple  and  non-simple  points  alternate  within  the  strip. 
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1 t 


Ap=l . 5h 


Figure  2.9.  A choice  of  A p = 1.5 h > h leads  to  strips  of  orientation  = 0 having 
simple  points  throughout  their  length. 


os  (0) 


Figure  2.10.  If  A p = 1.5 h > 2/icos(7t/4)  is  chosen,  strips  representing  three  DSLs 
with  angle  x/4  appear.  These  strips  alternate  with  strips  representing  two  DSLs  with 
the  same  angle. 


34 


* i 

^ ^Ap=0 . 5h 

Figure  2.11.  Illustration  of  the  case  where  A p = 0.5 h < h.  In  this  case  strips  having 
angle  9j  = 0 appear  with  no  points  belonging  to  them. 

For  Ap  < h max(|  cos  9\,  | sin  0|)  the  main  problem  is  that  the  digital  object 
defined  by  the  two  parallel  lines  may  not  be  8— connected  (see  Definition  1.4).  The 
fact  that  the  object  within  the  two  straight  lines  is  not  8— connected  means  that  there 
exist  columns  for  9 £ [0,  x/4)  (rows  for  9 £ [7t/4,  7t/2) ) that  have  no  pixel  belonging  to 
the  strip.  Another  problem  of  going  to  too  fine  of  a p— resolution,  is  that  there  will  be 
an  increase  in  the  of  spreading  of  the  DSL  segment  (see  Equation  2.3).  This  problem 
can  be  alleviated  by  postprocessing,  but  this  would  increase  the  computational  costs. 
Illustrative  examples  of  the  cases  where  A p < h max(|  cos  6\,  | sin  0|)  arise  when  9 — 0 
or  7t/4.  For  the  case  when  9 = 0 if  A p < h max(|  cos  9\,  | sin  0|)  = h there  will  be 
strips  that  have  no  pixels  at  all  (see  Figure  2.11).  The  same  applies  for  the  case 
9 = 7t/4  and  A p < \/2h/2. 

Summarizing,  from  the  perspective  of  Digital  Topology  theory  the  optimum 
solution  for  the  choice  of  Ap  is  an  adaptive  step  equal  to  Ap  = h max(|  cos  0|,  |sin0|) 
where  (j)  is  the  angle  of  the  strip.  For  the  cases  where  9 is  chosen  to  be  constant, 
there  would  be  objects  described  by  (pi,9j)  cells  that  are  either  not  8-connected 
(when  Ap  < h max(|  cos  0|,  | sin  0|))  or  have  simple  points  (when  A p > h max(|  cos  9 1,- 
I sin  <9|)) 

2.6  Performance  Prediction 


The  performance  prediction  problem  can  be  stated  as  follows  [60]: 
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Problem  2.1  Given  that  the  noise  has  a uniform  probability  in  the  image  space,  de- 
termine the  probability  density  function  in  the  parameter  space. 

Maitre  in  his  work  studied  the  effects  of  uniform  noise  for  both  ( p , 9)  and  (c,  m) 
parameterizations  and  both  circular  and  rectangular  images.  He  assumed  continuous 
image  and  transform  spaces  and  derived  the  probability  distributions  of  the  parameter 
spaces.  His  work  reached  conclusions  similar  to  those  of  an  earlier  work  published 
by  Alagar  and  Thiel  [4].  Hu  and  Destine  studied  the  same  effects  on  continuous 
image  and  parameter  space  [32].  Costa  and  Sandler  proposed  a nonuniform  A p 
distribution  to  reduce  the  nonuniformity  of  the  effects  of  the  noise  in  the  parameter- 
space  [15].  Cohen  and  Toussaint  [13]  and  Hu  and  Ma  [34]  proposed  postprocessing 
and  preprocessing  methods  to  reduce  the  effects  of  these  nonuniformities.  In  contrast 
to  all  previous  studies  we  will  be  considering  the  more  applicable  case  of  discrete 
image  and  parameter  spaces.  We  will  consider  the  line  parameterization  of  Duda  and 
Hart  (see  Equation  1.3)  and  we  will  assume  a rectangular  image  (see  Figure  1.2). 
The  point  of  origin  coincides  with  the  center  of  the  image. 

The  spatial  uniformity  of  the  noise  leads  us  to  the  conclusion  that  each  point 
has  a p probability  of  being  a noisy  point  and  q = 1 — p probability  of  not  being 
a noisy  point.  This  is  a Bernoulli  distribution  [29].  If  we  sample  r points  from  the 
image  space  (r  is  the  number  of  pixels  belonging  to  a strip)  it  is  easy  to  show  that 
the  probability  of  0 < n < r of  them  to  be  noisy  follows  a binomial  distribution  [29, 
page  103]: 

B(n,r,p)  = (^)p"(1  — PY~n-  (2.8) 

We  now  proceed  to  calculate  r for  a given  (p,  9)  for  a given  parameterization.  We  will 
use  simple  geometry  to  find  the  area  an  arbitrary  strip  (member  of  the  parameter 
space)  has  in  common  with  the  rectangular  image.  Obviously  there  is  a symmetry 
for  the  four  quadrants  of  the  image  so  we  only  need  to  consider  the  cases  where 
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Pi  E [0, 7V/2]  and  9 E (0, 7r/2)  (see  Figure  2.12).  A straight  line  with  parameters 
within  this  range  intersects  the  bounding  lines  of  the  image  at  the  points  (see  Figure 
2.13): 

• for  x = — N/2,  A{— N/2,  (pi  + N cos  9/2)/ sin  9), 

• for  x = N/2,  C(N/2,  (pi  — N cos  9/2)/ sin  9), 

• for  y = — N/2,  D((pi  + N sin  9/2)/  cos  9,  — N/2), 


• for  y = N/2,  B((pi  — N sin  9/2)/ cos  9,  N/2). 


Since  9 E (0,7t/2)  this  is  a decreasing  function  (<f>  E (ff/2, 7r)).  Therefore,  the  cases 
that  AB  or  CD  belong  to  the  image  are  excluded.  Hence,  the  segment  of  the  straight 
line  that  belongs  to  the  image  starts  with  either  point  A or  B and  ends  with  either 
point  C or  D.  In  order  for  B to  belong  to  the  boundary  of  the  image  the  following 
should  hold: 


-N  pi- %sm  9 N 

< — 2 < __ 

2 cos  9 2 


(2.9) 


Solving  for  pi  we  get 


2"(sin  9 


N 

cos  9)  < pi  < —(cos  9 + sin  9). 


(2.10) 


If  we  do  the  same  for  point  C we  get  that  in  order  for  that  point  to  belong  to  the 
boundary  of  the  image  the  following  should  hold: 


for  point  A, 


and  for  point  D, 


N . N 

—(cos  9 — sin  9)  < pi  < —(cos  9 + sin  9 ), 

2 


N . N 

— (—  cos  9 — sin  9)  < p,  < —(sin  9 — cos  9), 


N . N 

— (—  cos  9 — sin#)  < pi  < —(cos  9 — sin 9). 
z z 


(2.11) 


(2.12) 


(2.13) 
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For  9 £ (0,  tt/2)  and  0 < p,  < N/2  we  get  the  following  three  cases: 

1.  for  pi  > N/2\  cos  6 — sin0|  the  segment  BC  belongs  to  the  image, 

2.  for  N / 2(cos  9 — sin  9)  > p,-  > 7V/2(sin  9 — cos  9)  the  segment  BD  belongs  to  the 
image, 

3.  for  7V/2(sin  9 — cos  9)  > pi  > N /2(cos  9 — sin  9)  the  segment  AC  belongs  to  the 
image. 

The  lengths  of  the  segments  are  \BC\  = |pt  - N(cos  9 + sin  6>)/2|/|  sin0cos0|,  \BD\  = 
N/  cos  9 and  \AC\  = N/  sin0.  The  number  of  pixels  r belonging  to  a strip  is  given 
by: 

r{pi,9 ) = l(pi,9)w(pi,9)(lp™2el ),  (2.14) 

where  l(pi,9)  is  the  length  of  the  strip  inside  the  image  and  w(pi,9)  is  the  width  of 
that  strip.  If  we  consider  the  previous  cases  and  substitute  the  l(pi,9 ) we  calculated 
and  w(pi,  9)  --  A p = h we  will  get 

' k‘Ty.(in“^;n*)/i|  forft>iV|coS9-Sm«|/2, 

r=<  N/cos9  for  N(cos9  - s\n9)/2  > pi  > N(sin9  - cos0)/2,  , , 

N/  sin  9 for  A(sin  9 — cos  9)/2  > Pi  > N(cos  9 — sin  9) /2,  ^'15^ 

_ N for  9 = 0 or  9 = 7r/2. 

Hence,  the  number  of  votes  for  each  member  of  the  parameter  space  for  a uniformly 
corrupted  image  with  p probability  of  “success”  follows  a binomial  distribution  with 
r number  of  trials  (pixels)  where  r is  given  by  Equation  2.15  and  p percentage  of 
“success”  (success  means  that  a pixel  is  noisy). 

A plot  of  the  expected  values  of  the  accumulation  matrix  C for  a uniform  spatial 
distribution  of  the  noisy  pixels  in  the  image  domain  is  presented  in  Figure  2.14.  These 
values  were  calculated  using  the  mean  value  for  a binomial  distribution 


B — rp 


(2.16) 
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Figure  2.12.  Demonstration  of  the  symmetry  of  the  straight  lines  around  the  point 
of  origin. 


Figure  2.13.  Points  of  intersection  of  a line  with  0 < 9 < 7t/2  and  N/2  > p;  > 0 with 
the  straight  lines  that  form  the  boundaries  of  the  square  image. 
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Figure  2.14.  A plot  of  the  Hough  space  in  case  of  uniform  random  noise  for  a square 
image.  The  plot  represents  the  expected  values  for  the  first  quadrant:  6j  G [0,7t/2] 
and  pi  E [0,  N/2].  The  size  of  the  image  is  256  by  256  pixels. 

with  the  r calculated  from  Equation  2.15  and  p = 1.  The  value  of  p was  1 (all  pixels 
were  noisy). 

2.6.1  Experimental  Results 

In  Appendix  B we  provide  details  about  the  statistical  tests  we  performed  to 
show  that  the  votes  C(i,j)  follow  a binomial  distribution  for  a uniform  noise  distri- 
bution in  the  image  space. 

Plots  of  the  expected  values  according  to  the  theory  presented  (Equations  2.15 
and  2.16)  for  the  case  of  p = 1 uniform  noise  superimposed  on  the  experimentally 
obtained  values  are  presented  in  Figure  2.15.  Figures  2.15(a)  and  2.15(b)  show  two 
typical  cases.  For  the  first  case,  for  an  arbitrary  angle,  the  theoretical  and  the 
experimental  values  closely  match.  For  the  second  case,  for  6 = 7t/4,  due  to  the 
digitization  large  discrepancies  occur.  We  can  see  that,  in  the  mean,  the  predictions 
are  correct,  but  due  to  the  digitization  errors  we  get  these  large  discrepancies.  Figure 
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2.15(c)  which  shows  a vote  profile  for  the  p— direction  gives  a clear  picture  of  what 
is  going  on.  It  shows  that  for  most  of  the  strips  the  theory  works  well  but  there  are 
a number  of  orientations  for  which  strong  differences  appear. 

We  present  two  examples  to  demonstrate  the  correctness  of  the  theoretical  part 
of  this  section.  In  Figure  2.16(a)  we  show  a synthetic  image  corrupted  with  10% 
spatially  uniform  noise.  The  size  of  the  image  is  256  by  256  pixels.  Figure  2.16(b) 
shows  the  reconstructed  image  after  Hough  transform,  filtering  and  peak  detection 
and  reconstruction  of  the  detected  straight  lines.  The  detected  peaks  were  found  to 
have  angles  equal  to  7r/4  or  37t/4.  The  explanation  for  this  is  easy  if  one  considers 
what  is  shown  in  Figure  2.15(b).  For  the  strips  of  angles  7t/4  and  37r/4  large  fluc- 
tuations of  the  number  of  pixels  occur.  Due  to  the  spatial  uniformity  of  noise,  the 
number  of  votes  a strip  takes  is  proportional  to  the  number  of  pixels  belonging  to 
this  strip.  Therefore  the  large  fluctuations  of  the  number  of  pixels  belonging  to  these 
strips  are  the  reason  for  the  creation  of  local  maxima  at  cells  of  angles  7r/4  and  37t/4. 

Figure  2.17(a)  shows  the  image  of  a house.  The  size  of  the  image  is  512  by  512 
pixels.  In  Figure  2.17(b)  we  show  the  edge  detected  image,  while  in  Figure  2.17(c) 
we  show  the  detected  lines  superimposed  on  the  original  image.  Obviously  the  line 
pointed  to  by  the  arrow  does  not  represent  any  significant  linear  feature  and  since  its 
angle  is  equal  to  /4  it  is  attributed  to  the  phenomenon  discussed  previously.  The 
same  applies  for  the  rest  of  the  lines  shown  in  the  same  figure  and  are  parallel  to  the 
one  pointed  to  by  the  arrow.  Similar  situations  appeared  often  when  the  processed 
image  contained  large  concentrations  of  feature  points  (i.e.  texture)  in  some  part  of 
it.  In  the  image  we  are  discussing  the  cause  of  these  lines  is  the  leaves  of  the  trees 
that  appear  at  the  upper  right  part  of  the  image.  Suggested  solutions  to  this  problem 
are  the  following: 


1.  ignore  peaks  appearing  at  these  two  orientations 
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(a) 


(b) 


(c) 

Figure  2.15.  Three  profiles  comparing  the  performance  prediction  of  Equation  2.15 
with  the  results  obtained  with  the  use  of  the  Hough  transform,  (a)  For  9j  = 89;  (b) 
for  Qj  = 65;  and  (c)  for  pi  = 73.  For  Oj  = 65  (arrow  (i))  there  is  a large  discrepancy 
between  the  predicted  and  actual  value  due  to  the  digitization  errors.  The  same 
effect,  but  to  a lesser  extent,  appears  for  the  angles  pointed  to  by  arrow  (ii). 
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(a)  (b) 


Figure  2.16.  An  example  with  synthetic  images  that  demonstrates  the  danger  of 
recognizing  non-existing  lines  in  highly  noisy  environments,  (a)  A 256  by  256  syn- 
thetic image  corrupted  with  10%  uniform  noise;  (b)  lines  detected  after  the  Hough 
transform.  All  the  detected  lines  had  orientations  of  7r/4  and  37t/4. 

2.  choose  a quantization  step  that  excludes  these  two  orientations, 

3.  use  a ^-quantization  containing  the  two  orientations  but  start  the  quantization 
from  A0/2,  i.e.,  the  members  of  the  quantization  set  are:  A#/2,  3A0/2,  5A0/2, 
...,  7r  — AQ/2.  Using  this  strategy  we  avoid  having  7t/4  and  37t/4  in  the 
quantization  set, 

4.  use  the  quantization  proposed  in  Section  2.5. 

We  need  to  notice  that  this  phenomenon  does  not  only  appear  for  the  cases  of  the 
two  orientations,  7t/4  and  37t/4,  but  for  other  cases  as  well.  One  can  verify  this  by 
considering  the  strips  (ii)  of  Figure  2.15(c).  The  large  fluctuations  in  the  number  of 
pixels  they  have  may  give  highly  voted  local  maxima,  but  if  one  considers  the  strength 
of  these  fluctuations  compared  to  the  strength  of  the  fluctuations  of  the  strips  with 
angles  7t/4  and  37r/4  the  danger  is  significantly  smaller. 
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(a) 

Figure  2.17.  Real  image  demonstrating  the  danger  of  detecting  non-existing  lines, 
(a)  Image  “house”  obtained  from  INRIA’s  image  database;  (b)  detected  edges;  (c) 
detected  straight  lines  lines  superimposed  on  the  original  image; 


44 


0>) 

Figure  2.17.  continued 
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(c) 

Figure  2.17.  continued 


CHAPTER  3 

LENGTH  CALCULATION  USING  THE  HOUGH  TRANSFORM 


3.1  Introduction 

It  has  been  recognized  that  among  the  main  disadvantages  of  the  Hough  trans- 
form is  its  failure  to  provide  information  about  the  position  and  the  length  of  the 
digital  straight  line  segment  [3].  In  this  chapter  an  algorithm  is  proposed  for  the 
accurate  calculation  of  the  length  of  a digital  segment  by  its  Hough  transform.  Sim- 
ulations prove  the  superiority  of  the  new  technique  over  the  one  proposed  by  Akhtar 
and  Attiquzzaman  [3].  We  also  extend  the  proposed  technique  to  cases  where  noise 
and  other  objects  are  present  in  the  image.  This  represents  a significant  improvement 
over  recently  proposed  techniques  [3,  6]  which  relied  on  the  unrealistic  assumptions 
that 

• the  digital  straight  line  under  consideration  is  perfect, 

• there  are  no  other  objects  in  the  image,  and 

• no  noise  is  present. 

We  repeat,  here,  that  when  we  are  referring  to  an  “accurate”  length  estimator  we 
imply  the  distance  between  the  two  endpoints  of  the  DSL  segment.  This  estimator, 
referred  to  as  the  Most  Probable  Original  estimator,  is  not  optimal,  but  gives  results 
very  close  to  the  optimal  Best  Linear  Unbiased  estimator  [22].  We  also  have  to  add 
that  this  estimator  was  referred  to  as  “accurate”  length  in  an  earlier  work  [40]. 

3.2  Previous  Work 

Recently,  Akhtar  and  Atiquzzaman  [3]  used  the  Hough  transform  to  determine 

the  length  of  a digital  straight  line  segment  on  an  image.  Their  method  was  based 
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on  the  number  of  cells  to  which  the  segment  is  spread  out.  They  used  van  Veen  and 
Groen’s  Equation  1.4  to  derive  the  formula: 

(nj  ~ l)Ap 


1 = 


sin(  2 + dj'PA0) 


(3.1) 


where  / is  the  length  of  the  segment,  n £ is  the  number  of  rows  the  spreading  out  of 
the  votes  extends  to  for  the  j- th  segmentation  angle  (column  j of  the  accumulation 
matrix)  and  djiP  = | j — p|,  where  p — 8p/ A8  is  the  column  of  the  cell  in  which  the 
peak  is  found.  They  estimated  the  length  / using  several  angles  (columns  of  the 
accumulation  matrix)  and  they  averaged  these  estimations  to  improve  the  accuracy. 


3.3  A New  Algorithm  for  Length  Estimation 


3.3.1  Theory  and  Algorithm 


If  for  a DSL  segment  there  is  no  spreading,  the  votes  of  the  strip  where  the 
maximum  occurs  would  be  equal  to  the  number  of  pixels  that  belong  to  the  DSL 
segment.  The  spreading  of  the  DSL  segment  to  several  strips  is  almost  unavoidable 
so  the  votes  of  the  digital  segment  are  shared  among  several  cells  that  have  the  same 
orientation  (belonging  to  the  same  column  of  the  accumulation  matrix)  (see  Figure 
1.3).  If  vi  is  the  number  of  pixels  that  belong  to  the  digital  segment,  the  following 
relation  should  hold: 

VL=i2C(i,j)  (3.2) 

;=o 

where  N is  the  number  of  rows  of  the  vote  accumulation  matrix.  In  practice  it  is 
unecessary  to  sum  over  the  whole  interval  since  the  digital  segment  spreads  its  votes 
to  only  a small  number  of  cells.  Therefore  Equation  3.2  can  be  replaced  by 

lmax 

vl=  £ C(i,j)  (3.3) 

*=®min 

where  im;n  and  zmax  are  the  rows  of  the  extremal  strips  that  take  votes  from  the 
digital  segment  for  the  jf-th  column.  Relations  3.2  and  3.3  hold  for  all  j’s.  The  only 
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difference  is  that,  for  Equation  3.3,  as  we  go  away  from  the  column  of  the  primary  set 
the  summation  extends  to  more  and  more  members  of  the  column  (see  Figure  3.1). 
This  effect  happens  because  the  angle  between  the  digital  segment  and  the  strip  that 
is  represented  by  the  cell  (i,j)  becomes  larger  and  this  results  in  the  number  of  votes 
each  strip  takes  becoming  smaller.  In  Appendix  B we  prove  that  the  length  of  a DSL 
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Figure  3.1.  Demonstration  of  the  increase  in  spreading  when  moving  away  from  the 
column  containing  the  peak. 


segment  is  given  by  the  following  equation: 


(vL  - 1 )h 


(3.4) 


max(|  cos  6 1,  | sin  0|) 

We  can  combine  Equations  3.3  and  3.4  to  calculate  the  length  of  the  digital  segment. 
Because  an  accurate  value  of  9 is  not  known,  we  approximate  it  with  the  angle  9P  of 
the  cell  where  the  peak  appears. 

We  are  ready  now  to  state  the  proposed  algorithm  for  the  calculation  of  the 
length  of  a digital  segment: 
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Algorithm  2 Length  Estimation,  Ideal  Case 

perform  the  Hough  transform  of  the  Image, 
find  the  peak  cell,  C(m,p), 

use  Equation  3.3  to  find  the  number  of  pixels  that  belong  to  the  digital  segment, 
use  Equation  3.4  with  6 = 9p  = pA6  to  estimate  the  length  of  the  digital  segment. 

end 


3.3.2  Error  Estimation 


In  this  subsection  we  will  derive  error  bounds  for  the  proposed  length  estimation 
technique.  We  will  assume  that  we  can  determine  the  number  of  the  votes  of  the 
DSL  segment  accurately.  Therefore,  the  only  uncertainty  comes  from  the  digitization 
errors  of  the  parameter  space.  For  9S  £ [0,  7t/4],  |cos#s|  > |sin#s|.  Hence,  Equation 
3.4  becomes  l = (v^  — l)h/ cos  9S,  which  gives 


dvL  (uL-l)sin  6sd6Sxh 
cos  9 s ' cos 2#s  ’ 

If  we  substitute  ( d9s)max  — A#/2  into  Equation  3.5  we  will  get 


dl 


(vL  - 1)  sin  9s^p- 
cos2  9 , 


h. 


If  we  substitute  / = ( vl  — 1 )h/  cos  9S  we  can  get 


dl  tan#sA# 

T = 2~ 


(3.5) 


(3-6) 


(3.7) 


The  quantity  dl/l  becomes  maximum  when  the  right  hand  side  of  the  equation  above 
becomes  maximum.  This  happens  when  9S  = 7r/4  which  gives  cos6fi  = sin#,..  After 
substituting  into  Equation  3.6  we  get 


(-)  = 
V ^ )max  — 


A 9_ 

~2~ 


(3.8) 
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The  same  value  of  8S  gives  the  maximum  for  dl: 


(dl) 


max 


y/2(vL  - 1)A 8h 

2 


(3-9) 


It  can  easily  be  shown  that  Equations  3.9  and  3.8  hold  for  DSL  segments  of  any 
orientation.  For  the  relatively  coarse  parameterization  of  A#  = 7r/64  for  a 256  by  256 
image,  Equation  3.8  gives  ( dl/l)max  « 0.025  while  for  Ad  = 7t/256,  ( dl/l)max  « 0.006. 

3.3.3  Experimental  Results 


The  algorithm  was  tested  and  compared  with  the  results  of  the  algorithm  pro- 
posed by  Akhtar  and  Atiquzzaman  [3].  All  images  were  256  by  256  in  size.  The 
p— resolution  was  set  to  A p — h.  For  the  algorithm  proposed  by  Akhtar  and  Attiquz- 
zaman  [3]  both  alternatives  were  implemented  (averaging  of  the  length  estimation 
using  four,  six,  eight  and  ten  columns  from  each  side  of  the  peak  and  averaging  of 
the  length  estimation  using  ten  columns  with  smaller  indices  from  the  index  of  the 
peak  cell),  and  the  best  estimate  is  reported.  The  images  used  for  the  experiments 
contained  a single  perfect  straight  line  segment  and  no  noise.  In  Table  3.1  results  are 
shown  for  a resolution  equal  to  7r / 64.  The  superiority  of  our  method  is  apparent. 
Table  3.2  shows  length  estimations  with  a 8— resolution  equal  to  7t/256.  The  per- 
formance of  the  algorithm  proposed  by  Akhtar  and  Atiquzzaman  improves  but  it  is 
still  inferior  to  the  performance  of  our  algorithm  with  a resolution  four  times  larger. 
As  A 8 becomes  smaller,  the  performance  of  the  proposed  method  improves  because 
the  estimation  of  the  angle  of  the  digital  straight  line  segment  is  more  accurate  (see 
Equation  3.8). 

It  is  worth  mentioning  that  the  algorithm  uses  the  peak  as  the  best  estimate  of 
the  parameters  of  the  DSL  segment.  As  explained  in  the  previous  chapter,  this  is  not 
always  the  case.  We  need  to  notice  though,  that  this  is  not  the  main  problem  of  the 
algorithm.  The  algorithm  would  do  a good  job  even  if  the  orienation  estimation  is  not 
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Table  3.1.  Comparison  of  the  performance  of  the  two  methods  for  calculating  the 
length  of  a digital  segment  for  A 8 = 7t/64  and  A p - h. 


Vi 

c2 

Akhtar 

Proposed  Method 

Exact 

(5,8) 

(167,123) 

175.0 

200.4 

198.7 

(109,33) 

(217,102) 

113.1 

128.7 

128.2 

(35,50) 

(160,109) 

122.9 

137.2 

138.2 

(48,92) 

(200,167) 

165.1 

167.1 

169.5 

(60,62) 

(88,80) 

28.6 

32.5 

33.3 

(7,20) 

(123,90) 

119.1 

134.1 

135.5 

that  accurate.  The  main  problem  of  the  algorithm  is  that  it  relies  on  the  unrealistic 
assumptions  mentioned  in  Section  3.1.  If  these  assumptions  do  not  hold  it  would  be 
impossible  to  find  the  strips  the  digital  segment  spreads  its  votes  to.  The  purpose  of 
the  following  section  is  to  extend  the  algorithm  to  realistic  cases. 

3.4  Extensions  to  Realistic  Cases 
3.4.1  General  Considerations 

The  length  estimation  technique  presented  in  the  previous  sections  relied  on  the 
unrealistic  assumption  that  the  only  object  present  in  the  image  is  the  DSL  object 
and  no  noise  is  present.  Moreover,  it  assumes  that  the  DSL  under  consideration  is 
a perfect  one  according  to  Digital  Topology  theory  (see  Definition  1.8  and  Lemma 


Table  3.2.  Comparison  of  the  performance  of  the  two  methods  for  calculating  the 
length  of  a digital  segment  for  A 9 = 7t/256  and  A p = h. 


Vi 

c2 

Akhtar 

Proposed  Method 

Exact 

(5,8) 

(167,123) 

198.8 

196.9 

198.7 

(109,33) 

(217,102) 

124.5 

126.6 

128.2 

(35,50) 

(160,109) 

135.5 

137.2 

138.2 

(48,92) 

(200,167) 

153.3 

169.1 

169.5 

(60,62) 

(88,80) 

40.1 

32.7 

33.3 

(7,20) 

(123,90) 

127.1 

134.1 

135.5 
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1.1).  Of  course  this  case  never  appears  in  applications  and,  therefore,  there  is  a need 
for  an  extension  of  the  technique  to  nonideal  cases  where 

1.  noise  is  present, 

2.  other  objects  are  present,  and 

3.  the  edges  may  be  broken  and  noisy  in  the  sense  that  edge  pixels  do  not  form 
perfect  DSL  segments  according  to  Digital  Topology. 

In  this  section  we  will  extend  the  results  of  the  previous  section  to  the  cases 
where  the  image  is  corrupted  with  noise  and  other  objects  are  present.  Discussion 
on  the  extension  of  the  method  to  take  into  account  the  third  case  will  appear  at 
the  end  of  the  section  and  some  more  treatment  of  this  difficulty  will  appear  in  the 
chapter  that  follows.  We  present  tests  of  the  technique  with  synthetic  images:  digital 
straight  lines  built  with  the  use  of  the  Silicon  Graphics/Graphics  Library  or  the  use 
of  Bressenham’s  algorithm  [26,  57]  coded  by  the  author.  Tests  with  real  data  are 
presented  in  Chapter  4. 

3.4.2  Theory 

Suppose  that  a strip  (i,  j)  represents  best  the  parameters  of  a digital  straight 
line  segment.  We  will  be  referring  to  this  strip  as  “the  peak”  even  though  we  have 
not  shown  a method  to  determine  it.  The  votes  C{i,j')  this  strip  gets  come  from 
pixels  of  the  digital  segment  and  from  pixels  that  belong  to  other  objects  or  noise 
(see  Figure  3.2).  Therefore,  the  number  of  votes  of  the  peak  is  equal  to 

C(i,j)  = vL(i,j)  + vn(i,j)  (3.10) 

where  vl  is  the  number  of  votes  due  to  the  digital  straight  line  segment  and  vn 
is  the  number  of  votes  because  of  the  noise  and  other  objects.  Suppose  for  now 
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that  the  given  line  segment  spreads  its  votes  to  n^)  strips  and  the  number  of  these 
strips  is  odd.  We  can  also  suppose  that  these  strips  are  symmetrically  distributed 
around  the  peak.  Hence,  the  line  segment  spreads  its  votes  to  the  strips  C(i  — 


-,j)  and  by  summing  up  the  votes  of  these 


strips  we  get 


i+^i 


i+^ 


Y c(k,j)=  Y,  vL(k,j)+  Y vn(k,j). 


(3.11) 


k=i- 


;,4 

Solving  this  equation  for  £ \,j)  vL(k,j)  gives 

k=i-  np  \ -1 


»+-*- 


i+^ 


t+“ 


2 (fc.i)=  H C(k,j)~  Y ”n(k,j). 


(3.12) 


/c=t 


If  we  multiply  both  sides  of  this  equation  with  /*/  max(|  sin  Qj\,\  cos  0j\)  and  take  into 
account  Equation  3.4  we  obtain 


7 • 0 A . • 71 0 1 


k=i- 


k—i—  - 


max(  | sin  8j  | , | cos  6j  | ) 


-h, 


(3.13) 


where  fr’3'1  is  the  length  of  the  digital  straight  line  segment.  In  this  equation  there 
is  one  parameter  that  still  needs  to  be  evaluated:  the  votes  due  to  noise  and  other 
objects  vn(k,j).  For  the  estimation  of  votes  due  to  noise  and  other  objects  we  will  rely 
on  the  assumption  that  the  votes  due  to  objects  other  than  the  one  giving  the 
peak  are  distributed  among  the  primary  strips  and  their  neighboring  strips 
in  an  almost  uniform  manner.  A discussion  of  the  cases  when  this  assumption 
does  not  hold  is  presented  in  the  concluding  section  of  this  chapter. 

Since  we  “know”  that  the  digital  straight  line  segment  spreads  its  votes  to  the 
strips:  (k  — — -,j),  ...,(k+  np  '2  1,j),  the  strips  that  surround  them  are  voted 
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1 

Figure  3.2.  A figure  illustrating  the  various  sources  of  votes  for  a strip. 


only  by  noisy  pixels  and  pixels  from  other  objects.  Hence,  we  have 


n('u)  _ i n(y)  _ i 

1 ,j)  = vn(k--e— l J) 


(3.14) 


and 


n(*.i)  _ i n(*.i)  _ i 

C{k  H — — - 1-  1 ,j)  = vn(k  H — — 1-  l,i). 


(3.15) 


2 v 2 
So,  the  votes  of  these  two  strips  can  be  used  as  an  estimation  of  the  votes  due  to 
noisy  features  We  have 


n(hi)  _ i n(i,j)  _ 1 

C(k  — - 1 ,i)  + C(^  4 — — „ f l,j)  = 2 


(3.16) 


2 1 ’v”  1 2 
where  denotes  an  estimation  of  the  local  number  of  votes  due  to  noise  and  other 
objects,  local  meaning  for  neighboring  strips.  The  assumption  of  the  local  uniformity 
of  votes  due  to  noisy  features  gives: 


n(*.i)  _ 2 n(*.j)  _ i 

Vn  = Vn{k  1 ij)  = ...  = vn(k  H — h 1 ,j). 


(3.17) 
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Combining  Equations  3.13,  3.16  and  3.17  and  setting  /bJ)  _ we  get 

„£•»- 1 

1 ,4 

/(‘d)  = 


max( | cos  9j\,\  sin  #?•  | ) ^ ^ C(k,ri  1 

A;=,-?P-..  -1 

J 7,(*»d)  n(*>j)  — 1 1 

-£-(C(*  - -£-^—  - l,i)  + C(k  + + l,i)). 


(3.18) 


The  noise  estimation  can  be  made  more  accurate  if  we  use  two  more  strips  in  our 
noise  approximation  scheme  (one  from  each  side).  In  this  case  we  have 


/(m)  _ 


1 


i+^i 


max(|cos«,|,  |sin«,|)^  ^7,  1 

k=i-nP  „ * 


fil’d)  77,bd)  — t rr(®d)  — 1 

■— (C’(fc-^ U) + <?(*  + 2* 


+£(*- 


2 

nl‘d)  _ 1 


+ Ci)  + C(k  + 


2 

nl’d)  _ i 


~ 2,j)  + 

+ 2,i)). 


(3.19) 


Further  improvements  can  be  achieved  if  we  increase  the  number  of  strips  used  to 
estimate  the  noise.  In  this  case,  however,  we  run  the  risk  of  violating  the  assumption 
of  the  uniformity  of  votes  due  to  noisy  features. 

3.4.3  Algorithm 


The  algorithm  starts  with  the  assumption  that  every  member  of  the  parame- 
terization set  represents  a DSL  segment.  Hence,  Equation  3.18  or  3.19  is  used  to 
calculate  its  length.  It  is  easy  to  see  that  the  operation  of  using  Equation  3.18  to 
calculate  the  length  of  the  DSL  segment  corresponds  to  the  following  operation: 


— t 0 C T K 


(3.20) 
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where  the  symbol  ® denotes  filtering.  The  filter  t is  given  by: 


(*.j)  i 


t(ij)  - 


max( | cos  Qj\,\  sin  Oj | ) 


(3.21) 


where  the  filter  has  a length  equal  to  rajW)  + 2 and  is  to  be  slid  along  columns.  The 
arbitrary  member  of  matrix  K is  given  by: 

1 


K(i,j)  = — 


(3.22) 


max(  | sin  9 3\,\  cos  93  | ) 

The  operation  of  using  Equation  3.19  to  calculate  the  length  of  the  DSL  segments 
corresponds  to  Equation  3.20  where  t is  given  by: 

(i,j) 


t (,J) 


max(|  cos  9j |,  | sin  ^ | ) 


„(&) 

IP 


J.i) 


(3.23) 


In  order  to  apply  these  Liters  an  initial  approximation  of  njj,' J)  is  needed.  We  will 
assume  that  the  segment  has  the  maximum  possible  length  /£’J)  = (ln)max  = \/2 N 
for  an  N by  N image.  We  can  substitute  = l^’^h  into  Equation  2.3  to  calculate 
n{p'3)-  If  n is  even,  we  make  it  odd  so  that  the  Liter  is  symmetric  around  the 
central  cell.  An  alternative  approach,  when  n ^ is  even,  would  be  to  Liter  the  image 
with  two  operators  (see  Figure  3.3)  and  choose  the  largest  value.  After  the  Lltering,  a 
new  estimate  of  is  obtained.  If  we  substitute  the  new  /bd)  — Jnto  Equation 

2.3  we  will  Lnd  a new  nbJL  This  process  can  be  used  iteratively  until  no  change 
of  /(«)  occurs.  At  the  end  of  this  iterative  process  we  can  proceed  and  threshold 
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the  parameter  space.  The  steps  of  the  method  are  illustrated  in  the  algorithm  that 

follows: 

Algorithm  3 Length  Estimation,  Noisy  Case 

( ln)max  = V%N  /*  TV  by  N image  */ 

calculate  np  with  ln  = (ln)max  (see  Equation  2.3) 

np  — nP  + (nP  + 1)  mod  2 /*  make  filter  of  odd  length  */ 

derive  t from  Equation  3.21  or  3.23 

FC  = t ® C + K 

FC0id  — 0 
while  FC  ^ FC0id 
FCold  = FC 

calculate  np  (see  Equation  2.3) 

np  — np  + (np  + 1)  mod  2 /*  make  filter  of  odd  length  */ 

derive  t from  Equation  3.21  or  3.23 
FC  = t © C + K 

end 

threshold(EC) 

end 

The  length  estimates  are  held  in  matrix  FC.  Matrix  FC0id  is  kept  for 
comparison.  For  the  calculation  of  the  coefficients  of  the  filter  t , Equations  3.21  or 
3.23  can  be  used.  The  member  ny^  of  matrix  np  holds  the  spreading  of  the  DSL 
segment  that  is  best  represented  by  the  strip  with  parameters  (i,j). 

The  main  problem  with  the  alternative  implementation,  presented  in  Figure 
3.3(b),  is  that  the  final  estimation  is  not  unbiased.  For  example  for  an  image  with 
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max 


(a)  (b) 

Figure  3.3.  A schematic  representation  of  the  two  filtering  techniques,  (a)  If  rA,J)  is 
even  it  is  made  odd  so  that  the  filter  is  symmetric  around  the  center  pixel;  (b)  the 
other  option  is  to  use  two  nonsymmetric  filters  and  choose  the  maximum  value. 


nothing  but  uniform  noise,  the  filtered  Hough  space  will  have  nonzero  mean  value. 
This  is  the  reason  why  from  now  on  we  will  only  be  dealing  with  the  first  filtering 
implementation,  i.e.  Figure  3.3(a).  More  discussion  on  the  performance  prediction 
of  the  filtered  space  is  presented  in  the  following  subsection. 

3.4.4  Performance  Prediction  for  the  Filtered  Space 

In  this  subsection  we  will  treat  the  problem  of  performance  prediction  for  the 
filtered  space,  for  the  filter  of  Equation  3.23.  Similar  equations  can  be  derived  for  the 
filter  of  Equation  3.21.  The  motivation  of  this  study  is  the  need  to  derive  probabilities 
of  detecting  non-existent  lines  and  missing  real  ones.  In  other  words  it  will  provide 
us  with  a greater  insight  on  the  selection  of  the  threshold.  Our  assumption  about 
spatially  uniform  noise  for  the  image  domain  enabled  us  to  show  that  the  votes,  n,  a 
cell  takes  due  to  noise  follows  a binomial  distribution  (see  Equation  2.8):  F?(n,r, p). 
The  probability  p of  a pixel  to  be  noisy  is  constant  if  we  assume  a spatially  uniform 
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distribution  of  noise.  The  population  number  r can  also  be  assumed  constant  locally 
(see  Equation  2.15  and  Figure  2.14).  Since  the  “votes”  of  the  filtered  space  are  linear 
combinations  of  the  votes  of  the  parameter  space, 


"P-1  | 2 

FC(i,j)=  Y akC(k,j)  + K(i,j)  (3.24) 


n£-l+2  

with  2np- 1 2 ak  — 0 and  the  mean  value  C(i,j)  = rp  is  locally  constant  it  is  easy 

to  show  that 


FC(i,j)  = 


(3.25) 


It  is  easy  also  to  show  that  for  the  standard  deviation  for  a member  of  the  filtered 
space,  for  the  filter  of  Equation,  3.23, 


aFC(i,j)  = 


1 


■ JnP  + nl/4yJrp(  1 - p). 


(3.26) 


max(  | cos  6j  \ , | sin  6j  | ) 

An  equivalent  expression  can  be  derived  for  the  filter  of  Equation  3.21.  The  prob- 
ability that  n is  the  number  of  votes  at  a cell  (i,j)  of  the  filtered  space  is  given 

by 


np+ 3 


SB(n)=  Y n%k....«o). 


(3.27) 


Enp+3 

*=0  “*"*=  = 


k=  0 


Given  that  all  the  events  P;(rq ) are  independent  of  each  other  and  that  P,  is  a binomial 
with  a given  population  r and  success  probability  p we  get 


SB(n)=  Y II  B(nk,r,p).  (3.28) 

EL'o+3  k=0 

A plot  of  SB  is  presented  in  Figure  3.4. 

This  type  of  analysis  can  be  performed  for  the  filtered  space  for  the  filter  pro- 
posed by  Princen  et  al.  [71].  On  the  other  hand  performance  prediction  for  the  filter 
proposed  by  Leavers  and  Boyce  [55]  is  more  involved  because  the  variables  that  are 
added  are  not  statistically  independent.  For  both  our  filter  and  the  filter  proposed 
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Figure  3.4.  A plot  of  the  theoretically  predicted  distribution  of  the  filtered  Hough 
space  for  the  case  of  a uniform  noise  corrupted  image.  The  percentage  of  the  noise 
is  1%.  The  image  size  is  256  by  256  pixels.  The  strip  for  which  the  distribution  is 
plotted  has  256  pixels. 

by  Princen  et  al.  the  number  of  votes  of  the  strips  that  are  added  are  statistically 
independent.  This  is  the  case  because,  if  So,  S\,. . . ,Sn-\  are  the  sets  of  pixels  be- 
longing to  the  n strips  whose  votes  are  added  we  can  easily  show,  since  these  strips 
are  parallel  and  nonoverlapping  that  for  any  i and  j with  0 < i,j  < n — 1: 

Si  n Sj  = 0 (3.29) 

For  the  filter  proposed  by  Leavers  and  Boyce  this  is  not  the  case. 

3.4.5  Error  Estimation 

In  this  subsection  we  will  provide  error  bounds  for  the  length  estimation  tech- 
nique. We  will  assume  again  that  the  detected  peak  provides  the  most  accurate 
estimates  of  the  parameters  of  the  DSL  segment.  We  will  also  assume  that  the  only 
error  of  the  number  of  votes  dv  comes  from  the  uncertainty  of  the  estimation  of  the 
votes  due  to  noise.  Again  we  will  consider  only  the  case  when  6S  G [0,7t/4]  where 
|cos0s|  > |sin0s|.  The  other  cases  are  symmetric  and  can  be  treated  in  the  same 
fashion.  We  have  v = (v^  — l ) + dv  where  v is  the  number  of  votes  of  the  peak  cell 
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of  the  filtered  space,  Vi  the  number  of  votes  of  the  DSL  segment  (pixels),  and  dv  the 
error  in  the  estimation  of  the  number  of  votes.  We  can  substitute  v into  Equation 
3.4  to  get  (after  expansion) 


dvL  (vL  - I)d9ssin9s 
dl  = (7T7W  + ZZiTZ )h- 


cos  0. 


cos2  9 „ 


(3.30) 


If  we  divide  the  above  by  ( vl  — 1 )h/  cos  9S  and  recognize  that  this  is  equal  to  the 
length  of  the  DSL  segment  we  get 

— = h tan  9sd9.  (3.31) 

/ vL-  1 


Taking  a closer  look  at  this  equation  we  can  see  that  the  first  term  accounts  for  errors 
due  to  the  noise  while  the  second  term  accounts  for  errors  due  to  the  discretization  of 
the  parameter  space.  We  can  see  that  as  d9  decreases  the  contribution  of  the  second 
term  decreases  linearly  to  become  0 for  9S  = 0.  The  largest  error  for  a given  length 
appears  for  9S  — 7t/4.  For  this  case 


dl, 

7 


dvL 
vl  - 1 


+ d9. 


(3.32) 


Figure  3.5(a)  shows  a plot  of  dl/l  versus  the  number  of  pixels  of  a straight  line  (one 
should  remember  that  as  the  number  of  pixels  of  the  straight  line  increases  so  does  the 
length  of  it).  The  quantity  dv  was  set  equal  to  a pc,  given  by  Equation  3.26,  while  d9 
was  set  equal  to  A9/2.  The  error  drops  as  the  length  of  the  DSL  segment  increases. 
The  bumps  that  appear  at  certain  points  are  due  to  the  increase  of  the  worst  case 
spreading  which  is  plotted  as  a function  of  the  length  in  Figure  3.5(b).  Even  though 
the  worst  case  spreading  does  not  appear  in  Equation  3.31,  dl/l  implicitly  depends 
on  it.  This  is  apparent  if  one  notices  that  dv  depends  on  the  worst  case  spreading 
(see  Equation  3.26). 
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(a) 


two  filters 

filter  of  odd  length 


(b) 

Figure  3.5.  Plot  of  the  error  as  a function  of  the  length  of  a DSL  segment.  The  noise 
level  is  assumed  to  be  equal  to  1%.  A 9 — x/256,  A p = h.  The  size  of  the  image 
was  256  by  256.  (a)  Error  as  a function  of  the  length.  The  continuous  line  shows  the 
error  as  estimated  in  the  presented  implementation.  The  dotted  line  shows  the  error 
presented  in  the  alternative  implementation  when  two  filters  are  used  for  the  cell  if 
the  worst  case  spreading  is  even,  (b)  Worst  case  spreading  as  a function  of  the 
length.  The  continuous  line  shows  the  real  value  of  the  worst  case  spreading,  while 
the  dotted  line  takes  into  account  the  fact  that  the  worst  case  spreading  is  increased 
by  one  when  it  is  even  (for  our  implementation  of  the  length  estimation  algorithm). 
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3.4.6  Experimental  Results 

Results  for  a number  of  synthetic  images  containing  a single  straight  line  segment 
and  uniform  noise  are  presented  in  Table  3.3.  One  can  notice  that  the  accuracy 
generally  tends  to  become  better  as  the  noise  level  decreases.  This  is  what  Equation 
3.31  predicts.  Also  we  notice  that  for  smaller  lengths  the  error  is  significantly  larger 
than  for  long  segments.  This  is  also  in  accordance  with  Equation  3.31. 


Table  3.3.  Length  Estimation  under  various  uniform  noise  levels  in  images  containing 
a single  DSL  segment.  The  chosen  parameterization  was:  A 9 = 7r / 256  and  A p = h. 


Vertices 

Length  Estimation 

Exact 

Vi 

v2 

1% 

2% 

3% 

4% 

5% 

10% 

(5,8) 

(167,123) 

190.94 

196.80 

200.81 

195.57 

195.64 

186.93 

198.7 

(109,33) 

(217,102) 

130.63 

127.74 

126.56 

129.14 

128.55 

125.86 

128.2 

(35,50) 

(160,109) 

142.15 

136.89 

135.79 

134.96 

134.13 

129.98 

138.2 

(48,92) 

(200,167) 

170.45 

174.09 

173.80 

179.69 

173.53 

168.03 

169.5 

(60,62) 

(88,80) 

40.46 

37.66 

37.97 

38.28 

40.15 

44.82 

33.3 

(7,20) 

(123,90) 

139.32 

140.14 

137.45 

138.92 

143.81 

132.86 

135.5 

3.4.7  Discussion 

Our  initial  goal  of  proposing  an  algorithm  that  would  determine  the  length  of 
DSL  segments  that  may  exist  in  a digital  image  proved  to  be  fruitful  in  more  ways 
than  one.  It  enabled  us  to  construct  a new  filter  with  which  to  convolve  the  parameter 
space,  and  to  determine  clusters  in  the  parameter  space  that  represent  DSL  segments 
in  the  image.  Compared  to  the  previously  proposed  filtering  techniques  the  new 
filtering  method  has  the  great  advantage  of  an  intuitive  selection  of  the  threshold. 
The  length  of  the  DSL  segments  that  one  wants  to  detect  is  the  threshold  for  the 
filtered  space.  It  is  easy  to  calculate  the  danger  of  detecting  non-existing  lines  can 
easily  be  calculated  if  one  considers  the  threshold,  the  noise  level  and  the  calculations 
provided  in  Subsection  3.4.4.  A disadvantage  of  the  method  is  that  the  same  amount 
of  noise  has  different  impact  in  different  directions,  i.e.  for  strips  with  angle  equal 
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to  7r/4  a noisy  pixel  will  have  more  impact  than  for  strips  with  angle  equal  to  0. 
This  is  the  case  because  of  the  dividing  factor  that  appears  in  the  right  hand  sides 
of  Equations  3.21  and  3.23.  The  price  we  pay  to  get  back  the  signal  (the  length  of 
the  DSL  segment)  and  to  set  the  mean  value  of  the  noise  to  zero,  is  aggravation  of 
the  noise  level  (see  Equation  3.26). 

The  nyrt  that  was  used  in  our  algorithm  is  an  overestimation  of  the  real  value. 
Extracting  the  real  value  of  the  spreading  of  a DSL  segment  can  be  done  only  in 
trivial  cases.  So,  we  have  no  option  but  to  rely  on  this  overestimation  which  in  effect 
increases  the  uncertainty  (standard  deviation)  of  estimation  of  the  length  of  the  DSL 
segment.  If  the  primitives  to  be  recognized  slightly  deviate  from  the  linear  model,  this 
deviation  will  be  severly  penalized  by  our  filtering  technique.  This  happens  because 
the  votes  of  highly  voted  strips  will  be  subtracted  after  being  multiplied  by  0(np). 
An  easy  way  to  solve  this  problem  is  to  make  the  length  of  the  window  of  the  strips 
to  be  added  larger  (i.e.  instead  of  np  make  it  np  + 2 ),  as  in  Figure  3.6.  It  should  be 
noted,  though,  that  this  strategy  will  lead  to  a filtered  space  with  higher  noise  level. 
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(b) 


Figure  3.6.  Illustration  of  the  method  that  makes  the  filtering  strategy  more  robust, 
(a)  Segment  Lq  is  spread  to  more  than  np  strips.  This  may  happen  because  it  is 
not  exactly  a linear  segment  or  because  of  suboptimal  edge  detection.  The  standard 
filtering  strategy  may  fail  to  detect  its  presence  because  of  the  negative  weight  it 
puts  to  strips  that  get  voted  by  it.  (b)  By  making  the  window  larger,  say  np  + 2 this 
problem  is  resolved. 


CHAPTER  4 

AN  ALGORITHM  FOR  CONVEX  POLYGON  RECOGNITION 


4.1  Introduction 

4.1.1  Recognition  of  Polygons 

The  recognition  of  polygonal  objects  with  the  use  of  a digital  computer  is  of 
paramount  importance  because  of  the  abundance  of  such  objects  in  industrial,  office 
and  even  outdoor  environments.  These  systems  can  be  applied  to  manufacturing 
[91],  vehicle  recognition  [44,  87],  inspection  [12,  79,  80]  and  building  recognition  from 
aerial  images  [1].  It  should  also  be  noticed  that,  often,  arbitrary  shapes  can  be 
approximated  by  polygons  [50,  81].  This  approximation  greatly  facilitates  the  task 
of  object  recognition  and  localization.  Therefore,  it  is  not  surprising  that  a great 
number  of  polygon  recognition  algorithms  have  been  developed. 

Ventura  and  Chen  [90]  used  a model  based  technique  to  match  nonconvex  poly- 
gons. Their  method  calculates  the  translation  and  orientation  of  the  polygon  by 
minimizing  the  sum  of  the  square  error  between  the  input  and  the  reference  data. 
In  a followup  paper  Ventura  et  al.  [91]  used  the  minimization  of  the  minimum  zone 
error  towards  the  matching  of  the  raw  data  with  the  reference  data.  Both  works 
were  more  oriented  towards  matching  and  less  at  recognizing  the  polygonal  objects. 
In  other  words  they  assumed  that  the  very  important  step  of  the  grouping  of  the 
primitives  had,  somehow,  been  solved  [89]. 

4.1.2  Grouping  of  Primitives 

We  review,  briefly,  the  grouping  primitives  that  are  of  importance  for  the  recog- 
nition scheme  that  we  will  be  presenting  in  the  following  section.  From  these,  con- 
vexity has  recently  been  recognized  as  an  important  means  towards  the  recognition 


66 


67 


of  objects  [37,  42],  Two  more  grouping  methods,  proximity  and  parallelism,  have 
long  been  regarded  as  very  important  for  object  recognition  both  by  humans  [8]  and 
by  the  computer  [58]. 

Convexity 

Grouping  primitives  (i.e.,  straight  line  segments)  with  the  use  of  convexity  has  re- 
ceived considerable  attention  in  the  past  few  years  [37,  42].  Huttenlocher  and  Wayner 
[37]  grouped  line  segment  using  as  criteria  the  proximity  between  their  vertices  and 
the  creation  of  convex  groups.  Jacobs’  approach  [42]  had  the  same  philosophy  (i.e., 
use  of  convexity  and  proximity)  for  a model-based  object  recognition  algorithm.  He 
showed  experimentally  and  theoretically  that  salient  convex  groups  are  very  unlikely 
to  appear  in  a scene  by  chance. 

Proximity 

The  role  of  proximity  for  the  task  of  object  recognition  is  of  great  importance  for 
humans  and  machines.  The  Gestalt  psychology  school,  that  had  great  influence  on 
the  areas  of  perception  and  perceptual  organization,  considers  proximity  as  one  of  the 
major  categories  for  grouping  of  features  [58].  According  to  the  proximity  grouping 
rule:  elements  that  are  closer  are  grouped  together.  The  law  of  proximity  has  found 
great  use  in  the  area  of  pattern  recognition  [37,  42,  58,  61].  It  is  very  unfortunate, 
though,  that  there  are  no  techniques  that  can  extract  proximity  information  from  the 
standard  Hough  transform.  Since  the  standard  Hough  transform  cannot  determine 
the  vertices  of  the  segments  it  detects,  it  cannot  provide  reliable  information  about 
their  proximity  [36]. 

Parallelism 

Parallelism  plays  an  important  role  for  the  object  recognition  by  humans  and  by 
computers  [8,  58].  Engelbrecht  and  Wahl  [25]  used  parallelism  towards  proposing 
an  industrial  object  recognition  scheme.  Parallelism  was  used  also  in  Lowe’s  object 
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recognition  system  [58].  Ioannou  and  Dugan  [41]  used  parallelism  in  the  Hough  space 
to  recognize  parallelograms  in  digital  images  under  noisy  conditions. 

4.1.3  Polygon  recognition  with  the  use  of  the  Hough  transform 

Engelbrecht  and  Wahl  [25]  were  the  first  to  lay  down  the  rules  for  the  recognition 
of  polyhedral  objects  with  the  use  of  the  Hough  transform.  Their  method  is  model 
based  and  used  a number  of  nonaccidental  properties,  such  as  parallelism,  collinearity, 
and  cotermination.  They  combined  these  properties  with  information  from  the  votes 
of  the  peaks  in  the  Hough  space  to  recognize  polyhedral  objects.  The  processing  took 
place  in  the  Hough  space.  Their  paper  clearly  states  a set  of  rules  a system  should 
use  towards  the  recognition  of  these  objects  but  does  not  present  any  computational 
implementation  of  such  a system.  Many  questions  still  remain  unanswered: 

1.  How  does  the  system  deal  with  noise ? Noise  was  shown  (see  Chapter  3)  to 
influence  the  number  of  votes  and,  therefore,  the  estimated  length  of  a DSL 
segment.  How  does  noise  affect  the  overall  performance  of  the  system?  How 
does  the  performance  of  the  system  changes  with  increasing  noise  level? 

2.  How  do  occlusions  affect  the  system? 

3.  How  does  the  system  deal  with  the  errors  due  to  the  digitization  of  the  image 
and  the  parameter  space?  It  has  been  shown  (see  Chapter  2)  that  due  to  the 
digitization  errors  the  exact  parameters  of  a DSL  segment  cannot  be  determined 
even  in  the  ideal  case.  Hence,  after  the  Hough  transform  we  end  up  with  noisy 
estimates  of  the  pararameters  of  a DSL  segment  (ps,ds).  It  is  easy  to  show 
that  these  errors  may  have  serious  effects  on  the  determination  of  corners  and 
junctions. 
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4.  It  has  been  recognized  that  junctions  are  difficult  to  recover  from  binary  images 
and  a great  deal  of  research  effort  has  been  done  towards  this  direction  [61]. 
No  computational  implementation  for  junction  detection  was  presented  in  this 
paper. 

A couple  of  methods  for  the  recognition  of  polygonal  shapes  have  been  presented 
by  Davies  [19,  pages  284-300].  The  first  method  uses  the  generalized  Hough  trans- 
form parameterization  of  the  straight  line  segment.  This  parameterization  has  the 
advantage  of  providing  information  about  the  positions  of  the  vertices  of  the  DSL 
segment  at  the  expense  of  tremendous  computational  and  storage  costs.  In  order  to 
reduce  the  costs  Davies  suggested  the  use  of  improved  Hough  transform  implementa- 
tions, such  as  the  “fast”  Hough  transform  [56]  and  the  “adaptive”  Hough  transform 
[38].  This,  as  he  points  out,  may  have  effects  on  the  performance  of  the  method, 
both  in  terms  of  detection  and  localization.  These  effects  are  aggravated  in  cluttered 
environments.  The  second  method  is  also  based  on  the  generalized  Hough  transform. 
The  scheme  can  detect  the  modeled  polygon  under  any  rotation  and  translation.  The 
algorithm  is  simple  and  noise  resistant  but  even  for  the  simple  case  (exact  shape  and 
orientation  of  polygon  known),  it  suffers  from  excessive  computational  and  storage 
requirements  [19]. 

4.2  An  Overview  of  our  Approach 

In  this  subsection  we  present  the  hypothesis  about  the  environment  we  are 
dealing  with,  and  the  requirements  of  our  system: 

1.  Task:  our  task  is  to  recognize  convex  polygons, 

2.  Environment: 


(a)  the  number  of  polygonal  objects  in  the  image  is  relatively  small  (i.e.,  1-4), 
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(b)  the  environment  may  be  complicated  in  the  sense  that  many  other  non- 
polygonal  objects  may  exist, 

(c)  occlusions  may  appear  but  they  do  not  occlude  a whole  side  of  a polygon, 

Our  solution  towards  this  task  will  be  the  use  of  a standard  Hough  transform- 
based  technique.  A schematic  diagram  of  our  algorithm  is  presented  in  Figure  4.1. 
The  advantages  of  our  approach  are  the  following: 

1.  Computational  and  storage  savings.  Compared  with  the  approach  proposed  by 
Davies  [19]  the  proposed  algorithm  is  extremely  efficient  in  terms  of  computa- 
tion and  storage. 

2.  Ease  of  implementation. 

3.  Occlusion  and  noise  resistance. 

4.  Data  driven  recognition.  Any  instance  of  the  convex  polygons  can  be  recog- 
nized. The  method  does  not  require  a model  of  the  polygon  to  be  recognized. 
Therefore,  it  is  possible  to  recognize  convex  polygons  under  any  translation, 
rotation  and/or  scaling. 

Drawbacks  of  the  method  are 

1.  Dependence  on  the  peak  detection.  If  a whole  side  of  the  polygon  is  occluded 
there  is  no  way  of  recovering  the  polygonal  object. 

2.  Computational  requirements.  If  the  environment  becomes  too  complicated,  in 
the  sense  that  many  straight  lines  exist  in  the  image,  the  computational  costs 
increase  rapidly. 


71 


Figure  4.1.  Schematic  diagram  of  the  proposed  algorithm. 
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As  we  have  already  said,  until  now  there  is  no  reliable  technique  that  can  extract 
proximity  information  from  the  Hough  space.  To  overcome  this  problem  we  combine 
convexity  and  length  information  for  the  recognition  of  convex  polygons. 

4.3  Theory  and  Implementation 

4.3.1  Edge  Detection 

Our  intention  was  to  keep  a system  as  independent  and  flexible  as  possible. 
There  are  various  application  tailored  edge  detection  techniques.  Some  edge  detec- 
tion methods  can  be  applied  optimally  towards  certain  tasks  but  they  may  be  too 
computational  expensive  when  a more  general  recognition  framework  is  required. 
Such  an  edge  detector  is  the  one  proposed  by  Perona  [68]  which  is  tailored  towards 
the  extraction  of  junctions.  His  method  is  based  on  a steerable-scalable  kernel  and 
requires  convolutions  with  32  filters. 

On  the  other  hand  accuracy  is  always  important  in  object  recognition  algo- 
rithms. This  is  also  the  case  for  our  system.  Inaccuracies  introduced  during  the  edge 
detection  will  have  two  negative  results: 

1.  Error  in  the  localization  of  the  object:  This  is  the  result  of  inaccuracies  intro- 
duced during  the  edge  detection  stage.  In  fact,  the  inaccuracies  in  the  object 
localization  are  increased  during  later  stages  of  the  algorithm  (segmentation 
errors  during  the  Hough  Transfrom). 

2.  Inaccurate  length  estimation:  Inaccuracies  of  the  edge  detection  may  increase 
the  spreading  of  DSL  segments.  This  obviously  has  effects  on  the  number  of 
votes  in  the  Hough  space. 

For  the  examples  with  real  data  presented  in  this  section,  the  Canny  edge  detector 
was  the  method  of  choice  [10].  The  code  that  we  used  was  written  by  Bob  Fisher, 
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Dave  Croft  and  Andrew  Fitzgibbon  of  the  Robot  and  Vision  Group,  Department  of 
Artificial  Intelligence,  University  of  Edinburgh.  The  edge  detection  was  followed  by 
thresholding  and  thinning  [84].  In  some  examples  we  used  the  Sobel  and  Roberts 
edge  detectors  provided  as  part  of  the  Image  Processing  Toolbox  of  MATLAB.  The 
effects  of  the  errors  due  to  the  edge  detection  on  the  accuracy  of  the  estimates  of  the 
Hough  transform  have  been  discussed  by  Palmer  et  al.  [65]  and  Matas  and  Kittler 
[61]. 

4.3.2  Hough  transform 

The  algorithm  described  in  the  first  chapter  was  implemented  as  a part  of  the 
proposed  system.  The  retina  we  used  was  circular.  Lines  lying  outside  the  retina  were 
not  taken  into  account  but  this  does  not  present  a problem  because  they  are  of  small 
length  and,  therefore,  not  important  (see  Figure  4.2).  The  algorithm  can  be  used 
with  a rectangular  retina  without  significant  modifications.  For  all  the  experiments 
the  size  of  the  parameter  space  was  chosen  to  be  the  same  as  the  image  space  (i.e.,  for 
a 512  by  512  image,  the  parameter  space  was  512  by  512).  This  implies  that  A p — h 
and  A 9 — t/N  where  A is  the  size  of  the  image.  The  center  of  the  image  was  chosen 
to  be  the  center  of  the  pixel  with  coordinates  (A/2,  A/2).  All  images  start  from 
(0,  0).  Hence,  the  image  is  slightly  nonsymmetric  with  respect  to  the  point  of  origin. 

4.3.3  Grouping  of  Primitives 

The  output  of  the  Hough  transform,  after  peak  detection,  is  a set  of  n straight 
line  segments  associated  with  their  estimated  length.  The  criterion  for  grouping  these 
primitives  is  the  creation  of  convex  groups.  Obviously  if  one  is  looking  for  triangles 
he  will  have  to  test  0(n3)  hypotheses  while  for  quadrilaterals  he  will  have  to  test 
0(n4)  hypotheses.  Any  set  of  three  straight  lines  defines  a triangle  and  also  any 
set  of  four  straight  lines  defines  a quadrilateral.  A more  detailed  discussion  of  this 
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A 

B 


Figure  4.2.  The  maximum  length  of  a straight  line  segment  lying  outside  the  circular 
retina  can  be  shown  to  be  equal  to  N/\/ 2 where  N by  N is  the  size  of  the  image. 

subject  is  presented  in  Chapter  5.  To  reduce  the  number  of  hypotheses  to  be  tested 
one  can  use  parallelism  and  size  (length  of  straight  line  segments)  information.  The 
use  of  parallelism  (if  applicable)  can  make  the  proposed  algorithm  very  efficient  by 
greatly  reducing  the  number  of  polygons  to  be  tested  (i.e.  parallel  segments  cannot 
belong  to  the  same  triangle,  parallel  segments  are  highly  likely  to  belong  to  the  same 
quadrilateral,  etc.).  Size  information  can  also  be  used  for  a reduction  of  the  search 
space.  An  example  is  that  if  the  lengths  of  three  segments  are  /i,  /2  and  /3,  from  the 
triangular  inequality: 

h — h < h-  (4-1) 

One  can  make  use  of  this  inequality  to  achieve  efficiency  in  the  grouping  process.  The 
problem  with  the  use  of  Equation  4.1  stems  from  the  fact  that  the  length  estimations 
are  not  accurate  at  all.  In  fact  the  large  discrepancies  that  may  occur  due  to  occlusions 
require  relaxing  of  the  above  rule. 

The  fact  that  the  computational  complexity  increases  combinatorially  with  in- 
creasing number  of  detected  line  segments  represents  a good  reason  for  an  accurate 
peak  localization.  During  the  peak  detection  stage  there  are  many  redundant  straight 
lines  (ghost  lines  [66]).  Chang  and  Hashimoto  reported  problems  with  such  lines  when 
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the  threshold  for  the  peak  detection  was  set  to  a too  low  value  [11].  If  for  each  peak 
there  are  two  straight  lines  (i.e.  one  ghost  line  per  line  segment)  the  computational 
costs  for  triangle  detection  would  increase  by  a factor  of  eight,  while  for  a quadri- 
lateral detection  they  would  increase  by  a factor  of  sixteen.  Efficient  grouping  was 
beyond  the  scope  of  this  dissertation  which  as  we  already  said  deals  with  relatively 
simple  environments. 

4.3.4  Matching  Coefficient 

If  one  relies  on  the  peak  detection  of  the  standard  Hough  transform  for  convex 
polygon  recognition  he  will  end  up  with  a large  number  of  “recognized”  polygons.  For 
example  if  n = 10  is  the  number  of  straight  lines  detected  by  the  Hough  transform, 
the  number  of  possible  triangles  that  may  exist  is  ( ^ j = 120.  We  therefore  need  a 
means  of  testing  the  existence  of  these  polygons.  This  section  discusses  the  use  of 
the  length  information  provided  by  the  Hough  transform  with  the  filtering  technique 
presented  in  the  previous  chapter  towards  the  task  of  verification  of  the  presence  of 
convex  polygons. 

Given  that  a set  of  n straight  lines  defines  a convex  n- gon,  we  prove  in  Chapter 
5 that  this  n-gon  is  uniquely  determined.  In  the  same  chapter  we  also  provide  an 
algorithm  that  constructs  the  n- gon  from  the  parameters  of  its  edges  in  O(nlogn) 
time.  From  the  vertices  of  the  n-gon  one  can  easily  calculate  the  lengths  of  its 
sides:  /q,  l9n_\.  On  the  other  hand,  the  Hough  transform  with  the  method  de- 

scribed in  the  previous  chapter  provides  a length  estimation  for  the  detected  segments: 
Iff , , . . . , A measure  of  discrepancy  between  the  geometry  based  estimations 

and  the  Hough  based  estimations  is  given  by  the  matching  coefficient  ( MC ): 

MC  = ^k=0il\~1^ 

En—  1 ig 
k=0  lk 


(4.2) 
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Under  the  assumption  that  the  parameters  of  the  straight  lines  are  exactly  deter- 
mined by  the  peaks  of  the  Hough  transform,  the  MC  gives  a measure  of  the  portion 
of  the  boundary  of  the  convex  polygon  that  is  visible.  Obviously  the  MC  varies  from 
0 to  -f  oo.  For  MC  = 0 we  have  a perfect  match  and  we  are  “sure”  that  the  object  is 
there.  As  MC  increases  our  confidence  on  the  existence  of  the  object  decreases.  We 
have  to  point  out  that  the  matching  coefficient  defined  by  Equation  4.2  resembles  the 
salience  ratio  used  by  Jacobs  in  his  convex  polygon  detection  scheme  [42]. 

4.4  Results 

4.4.1  Synthetic  Data 

Testing  the  algorithm  with  synthetic  data  is  important  because  it  allows  us  to 
reach  conclusions  about  the  accuracy  of  the  method.  It  also  gives  us  an  idea  on  how 
the  important  factors  (i.e.,  the  MC)  change  as  some  parameters  change.  We  are  very 
interested  to  see  how  the  matching  coefficient  changes  with  increasing  noise  level. 

Example  1 

Figure  4.3(a)  shows  a convex  4-gon  in  a 256  by  256  image.  The  image  is  cor- 
rupted with  1%  uniform  noise.  Figure  4.3(b)  shows  the  detected  straight  lines  and 
Figure  4.3(c)  shows  the  reconstructed  4-gon.  Figure  4.6  shows  a plot  of  the  match- 
ing coefficient  as  a function  of  the  noise  level.  As  expected,  the  matching  coefficient 
generally  tends  to  increase  with  the  noise  level. 

Example  2 

Figure  4.4(a)  shows  another  convex  4-gon  in  a 256  by  256  image.  The  image  is 
corrupted  with  1%  uniform  noise.  Figure  4.4(b)  shows  the  detected  straight  lines  and 
Figure  4.4(c)  shows  the  reconstructed  4-gon.  Figure  4.6  shows  a plot  of  the  matching 
coefficient  as  a function  of  the  noise  level.  The  same  conclusion  as  with  the  previous 
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(a) 


(b) 


(c) 

Figure  4.3.  Figure  illustrating  the  steps  of  the  algorithm  for  the  4-gon  of  Example  1. 
The  vertices  of  the  polygon  were  A(220, 190),  77(50,120),  (7(50,50)  and  77(180,50). 
(a)  Original  synthetic  image  (noise  level  1%);  (b)  detected  lines;  (c)  recognized  4-gon. 
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example  can  be  reached  that  as  the  noise  level  increases,  the  matching  coefficient 
generally  tends  to  increase.  We  also  see  that  the  matching  coefficient  is  larger  than 
the  previous  case.  The  reason  is  that  the  two  polygons  differ  in  size.  The  polygon  of 
this  example  is  larger  and  therefore  less  affected  by  noise  as  Equation  3.32  predicts. 

Example  3 

Figure  4.5(a)  shows  a convex  4-gon  with  50%  of  one  side  missing.  Such  experi- 
ments were  performed  in  order  to  study  the  effects  of  occlusions.  Figure  4.5(b)  shows 
the  detected  straight  lines  and  Figure  4.5(c)  shows  the  reconstructed  4-gon.  The 
MC  is  plotted  in  4.6.  We  can  see  that  due  to  “occlusions”  its  MC  is  much  larger 
than  in  the  previous  two  cases. 

Further  insight  on  the  advantages  and  drawbacks  of  the  proposed  scheme  will 
be  acquired  with  the  study  of  it  with  real  data. 

4.4.2  Real  Data 

Cylinder  on  a box 

The  image  of  4.7(a)  was  provided  by  Dr.  Maria  Petrou  of  the  University  of 
Surrey.  It  is  a 256  by  256  image.  Figure  4.7(b)  shows  the  edge  detected  image.  The 
edge  detection  method  used  was  the  Roberts  edge  detection  provided  as  a part  of  the 
Image  Processing  Toolbox  of  MATLAB.  Figure  4.7(c)  shows  the  recognized  straight 
lines  superimposed  on  the  original  image  and  Figure  4.7(d)  shows  the  recognized  4- 
gon  superimposed  on  the  original  image.  The  matching  coefficient  was  found  to  be 
equal  to  0.2237  The  second  lowest  matching  coefficient,  was  of  the  4-gon  of  Figure 
4.7(e).  Its  matching  coefficient  was  equal  to  0.2353. 

Figure  4.8(b)  shows  the  results  of  the  edge  detection  on  the  “cylinder  on  a box” 
image  for  the  same  edge  detector  but  for  a smaller  threshold  value.  We  performed  such 
experiments  in  order  to  show  how  the  system  responds  under  suboptimal  parameter 
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(a) 


(b) 


(c) 

Figure  4.4.  Figure  illustrating  the  steps  of  the  algorithm  for  the  4-gon  of  Example  2. 
The  vertices  of  the  4-gon  were  21(35,80),  5(80,100),  (7(130,110)  and  5(65,40).  (a) 
Original  synthetic  image  (noise  level  1%);  (b)  detected  lines;  (c)  recognized  4-gon. 
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(a) 


(b) 


(c) 

Figure  4.5.  Figure  illustrating  the  steps  of  the  algorithm  for  4-gon  of  Example  3. 
The  vertices  of  the  polygon  were  H(20, 10),  H(90,35),  (7(130,100)  and  T>(80, 110). 
Half  of  edge  connecting  vertices  A and  D was  missing,  (a)  Original  synthetic  image 
(noise  level  1%);  (b)  detected  lines;  (c)  recognized  4-gon. 
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Figure  4.6.  Matching  coefficients  for  the  4-gons  of  Examples  1-3. 

selection.  The  image  is  highly  noisy  especially  on  the  surface  of  the  cylinder.  Figure 
4.8(c)  shows  the  detected  lines  superimposed  on  the  original  image.  The  matching 
coefficient  for  the  real  4-gon  (see  Figure  4.8(d))  was  found  to  be  0.2615,  while  the 
matching  coefficient  for  three  ghost  4-gons  (see  Figure  4.8(e))  were  equal  to  0.1570, 
0.1584  and  0.2166.  The  fact  the  several  ghost  polygons  have  matching  coefficients 
less  than  the  matching  coefficient  of  the  real  polygon  is  attributed  to  the  combined 
effects  of  the  high  noise  levels  and  the  partial  occlusion  of  the  real  4-gon. 

Book  on  a table 

Figure  4.9(a)  shows  the  “Book  on  a table”  image.  This  is  a 512  by  512  image 
included  in  the  image  database  of  Carnegie  Mellon  University.  Figure  4.9(b)  shows 
the  results  of  the  Sobel  edge  detection.  Figure  4.9(c)  shows  the  straight  lines  that 
were  detected  through  the  use  of  the  Hough  transform.  The  recognized  4-gon  with 
the  lowest  MC  = 0.1391,  is  superimposed  on  the  original  image  in  Figure  4.9(d). 
Three  more  “ghost  “ 4-gons  are  superimposed  on  the  original  image  in  Figure  4.9(e) 
Their  matching  coefficients  were  0.1505,  0.1522  and  0.1895.  Had  the  right  edge  of 
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(c)  (d) 

Figure  4.7.  Figure  illustrating  the  steps  of  the  algorithm  on  the  image  “Cylinder  on 
a box”,  (a)  Original  image;  (b)  detected  edges;  (c)  detected  lines  superimposed  on 
the  original  image;  (d)  recognized  4-gon;  (e)  ghost  polygon. 
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(e) 


Figure  4.7.  continued 

the  real  4-gon,  especially  its  lower  part,  been  closer  to  a straight  line,  its  matching 
coefficient  would  have  been  much  smaller. 

Lab 

Figure  4.10(a)  shows  an  image  taken  in  the  Robotics  Laboratory  of  the  Depart- 
ment of  Nuclear  Engineering  Sciences  at  UF.  The  size  of  the  image  is  480  by  480. 
Figure  4.10(b)  shows  the  edge  detected  image  while  Figure  4.10(c)  shows  the  detected 
lines  superimposed  on  the  original  image.  Figure  4.10(d)  shows  two  detected  4-gons 
representing  the  same  real  one  (board  on  the  wall).  Their  matching  coefficients  were 
0.2201  and  0.2235.  The  ghost  parallelograms  shown  in  Figure  4.10(e)  had  matching 
coefficients  0.3120  and  0.3288. 

Figure  4.11(a)  shows  an  image  of  the  same  scene  with  the  only  difference  be- 
ing that  the  real  4-gon  is  partially  occluded  by  a plant.  Figure  4.11(b)  shows  the 
edge  detected  image  and  Figure  4.11(c)  shows  the  straight  lines  that  were  detected 
overlayed  on  the  original  image.  Two  ghost  4-gons  were  found  to  have  lower  MC s 
than  the  real  one  (Figure  4.11(d)).  The  matching  coefficients  for  the  ghost  4-gons 
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Figure  4.8.  Results  of  the  algorithm  on  the  image  “Cylinder  on  a box”  under  noisy 
conditions,  (a)  Original  image;  (b)  detected  edges;  (c)  detected  lines  superimposed 
on  the  original  image;  (d)  recognized  4-gon;  (e)  ghost  polygons. 
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(e) 

Figure  4.8.  continued 


(a) 

Figure  4.9.  Results  of  the  algorithm  on  the  image  “Book  on  a table”,  (a)  Original 
image;  (b)  detected  edges;  (c)  detected  lines  superimposed  on  the  original  image;  (d) 
recognized  4-gon;  (e)  ghost  polygons. 
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(c) 

Figure  4.9.  continued 
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(e) 

Figure  4.9.  continued 
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were  0.2593  and  0.2691.  while  the  matching  coefficient  for  the  real  4-gon  (see  Figure 
4.11(e))  was  equal  to  0.3041. 

Toys 

Figure  4.12(a)  shows  a 256  by  256  image  consisting  of  a number  of  polygons. 
This  image  was  provided  by  Dr.  Maria  Petrou  of  the  University  of  Surrey.  Obviously, 
our  system  would  have  great  difficulties  in  detecting  4-gons  in  such  an  environment 
because  of  accidental  matches.  The  large  number  of  straight  line  segments  would 
also  make  the  computational  requirements  extremely  high.  Therefore  the  task  for 
our  system  was  set  to  be  the  detection  of  the  hexagon  that  lies  in  the  middle  of  the 
image.  Figure  4.12(b)  shows  the  image  after  edge  detection  and  Figure  4.12(c)  shows 
the  straight  lines  that  were  detected  with  the  use  of  the  Hough  transform.  Figure 
4.12(d)  shows  the  recognized  hexagon.  Its  matching  coefficient  was  equal  to  0.1506. 
Figure  4.12(e)  shows  a ghost  hexagon  whose  matching  coeffient  was  equal  to  0.1596. 

4.5  Discussion 

The  applicability  of  our  length  estimation  technique  was  demonstrated  through 
its  use  in  a convex  polygon  recognition  scheme.  We  can  say  that  our  algorithm  is  not 
best  suited  in  applications  that  require  the  recognition  of  polygons  in  complicated 
environments  (large  number  of  polygons).  As  the  number  of  straight  lines  increases 
the  number  of  cases  to  be  explored  increases  combiatorially.  Also,  with  the  number 
of  lines  increasing  the  danger  of  accidental  matches  increases  too.  Therefore  for 
complicated  environments  more  sophisticated  recognition  techniques  would  be  more 
suitable  [27,  42],  On  the  other  hand  our  algorithm  performs  well  when  there  are  few 
polygons  to  be  recognized.  It  also  shows  robustness  to  noise  and  occlusions. 

The  choice  of  the  value  for  threshold  of  the  matching  coefficients  is  essential 
for  the  system.  If  too  low  a threshold  value  is  chosen  real  polygons  will  be  missed. 
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(a) 


Figure  4.10.  Recognition  of  the  board  on  the  image  “Lab”,  (a)  Original  Image;  (b) 
detected  edges;  (c)  detected  lines  superimposed  on  the  original  image;  (d)  recognized 
4-gon;  (e)  ghost  polygons. 
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(b) 


(c) 

Figure  4.10.  continued 


91 


(d) 


(e) 

Figure  4.10.  continued 
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(a) 


Figure  4.11.  Recognition  of  the  occluded  board  on  the  image  “Lab”,  (a)  Original 
image;  (b)  detected  edges;  (c)  detected  lines  superimposed  on  the  original  image;  (d) 
recognized  4-gon;  (e)  ghost  polygons. 
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(b) 


(c) 

Figure  4.11.  continued 
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(e) 

Figure  4.11.  continued 
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Figure  4.12.  Recognition  of  the  6-gon  on  the  image  “Toys”,  (a)  Original  image;  (b) 
detected  edges;  (c)  detected  lines  superimposed  on  the  original  image;  (d)  recognized 
6-gon;  (e)  ghost  hexagon. 
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(e) 


Figure  4.12.  continued 

On  the  other  hand  if  a large  value  is  chosen  we  will  end  up  with  a large  number  of 
ghost  polygons.  Empirical  results  suggest  that  for  unoccluded  4-gons  the  matching 
coefficient  rarely  becomes  larger  than  0.2.  Empirical  data,  again,  suggest  that  ghost 
polygons  with  a matching  coefficient  less  than  0.2  almost  never  appear.  Exceptions 
to  this  rule  were  found  for  the  cases  where  the  ghost  polygon  was  somehow  related 
to  a real  one  (see  Figure  4.9). 

These  findings  are  of  empirical  nature  and  cannot  be  stated  as  rules.  If  someone 
wants  to  form  a set  of  rules  concerning  the  choice  of  the  threshold  matching  coefficient 
he  has  to  consider  the  number  of  objects  in  the  scene,  the  size  of  the  objects,  possible 
occlusions,  etc. 


CHAPTER  5 

UNIQUENESS  OF  A CONVEX  POLYGON  IN  THE  HOUGH  DOMAIN 


5.1  Introduction 

In  this  chapter  we  present  a proof  of  the  uniqueness  of  the  representation  of  a 
convex  polygon  by  the  peaks  its  edges  give  in  the  Hough  domain.  We  will  assume 
an  idealized  model  (i.e.,  the  peaks  give  the  parameters  of  the  straight  lines  exactly). 
We  show,  through  a counter-example,  that  non-convex  polygons  are  not  uniquely 
determined  by  the  peaks  they  give  in  the  Hough  domain  (see  also  the  discussion  in 
[75]).  We  will  also  present  an  efficient  algorithm  that  reconstructs  a convex  polygon 
by  the  equations  of  its  edges. 

5.2  Background-Notation 

The  following  theorems  will  be  useful  as  we  develop  the  theory  and  the  algo- 
rithms in  this  chapter. 

Theorem  5.1  A ray  emanating  from  an  interior  point  of  a bounded  convex  figure  F 
intersects  the  boundary  of  F in  exactly  one  point. 

Proof:  See  [59,  page  6].  □ 

Theorem  5.2  A straight  line  cannot  have  more  than  two  common  points  with  the 
boundary  of  a convex  figure  F. 

Proof:  Easy  if  one  takes  into  account  Theorem  5.1.  □ 

Suppose  a straight  line  L and  a convex  figure  Q are  given.  There  are  five  possible 
arrangements  [59,  page  8]. 

• Case  1:  L and  Q have  no  points  in  common  (line  L0  in  Figure  5.1), 
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• Case  2:  L and  Q have  one  point  in  common  (line  L\  in  Figure  5.1), 

• Case  3:  L and  Q have  as  an  intersection  a line  segment  belonging  entirely  to 
the  boundary  of  Q (line  L2  in  Figure  5.1), 

• Case  4:  The  intersection  of  L and  Q lies  entirely  (with  the  exception  of  its 
endpoints)  inside  Q (line  L3  in  Figure  5.1), 

• Case  5:  (Only  for  unbounded  convex  objects)  The  intersection  of  L and  Q lies 
entirely  (with  the  exception  of  one  point)  inside  Q. 

We  state  two  properties  we  will  make  use  of,  for  the  proofs  of  Lemmas  5.1,  5.2, 
5.3  and  Theorem  5.3. 

Property  5. 1 a straight  line  cuts  the  plane  into  two  semi-planes  which  are  convex  sets. 

Property  5.2  if  the  intersection  of  two  convex  figures  is  not  the  empty  set,  it  is  convex. 

When  the  notation  f(A,B,...,P)  > 0 is  used,  it  is  meant  that  function  / becomes 
positive  for  all  points  A,  B , . . . , P.  Straight  lines  will  be  denoted  by  Li  and  their 
equations  by  ffix,  y)  — y — C{X  — mi  = 0.  The  notation  PqP\  . . . Pn-i  denotes  the  set 
of  vertices  of  a bounded  convex  n— gon  as  one  travels  its  boundary  counterclockwise 
(see  Figure  5.2).  The  notation  LqPqPi  . . . Pn_2Ln-i  denotes  an  unbounded  convex 
n— gon  as  one  travels  its  boundary  counterclockwise  (starts  from  ray  L0P0  and  ends 
with  ray  Pn_2Ln- 1,  see  Figure  5.3). 

5.3  The  uniqueness  of  the  Hough  transform  of  convex  polygons 
5.3.1  Main  Result 

In  order  to  prove  the  main  theorem  that  will  help  us  to  establish  the  uniqueness 
of  the  representation  of  convex  polygons  by  their  Hough  transformation,  Lemmas 
5.1,  5.2  and  5.3  will  be  needed. 
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Figure  5.1.  Illustration  of  the  four  different  relative  positions  of  a straight  line  with  a 
convex  shape:  line  Lq  has  no  common  points  with  the  convex  shape  Q,  line  L\  has  a 
single  common  point  with  the  boundary  of  Q , line  L 2 has  an  edge  segment  belonging 
to  the  boundary  of  Q and  line  L3  has  an  intersection  lying  entirely  inside  Q except 
for  its  endpoints. 


100 


* / 

/ 


Figure  5.2.  A bounded  convex  n— gon  which  is  cut  by  line  Ln.  Obviously  k vertices 
Po  . . . Ph-i  are  left  on  one  semi-plane  and  the  remaining  (n  — k ) vertices  of  the  n-gon 
are  on  the  other  semi-plane. 


Figure  5.3.  An  unbounded  convex  n— gon.  When  a straight  line  has  a common  point 
with  the  interior  of  the  n— gon  two  cases  exist:  case  (i)  line  Ln  cuts  the  unbounded 
n— gon  giving  one  bounded  polygon  and  one  unbounded  convex  polygon,  case  (ii) 
line  Ln  cuts  the  unbounded  n— gon  giving  two  unbounded  polygons. 
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Lemma  5.1  Suppose  P0Px...Pn_ 1 is  a bounded  convex  n-gon.  A straight  line  Ln,  which 
has  a common  point  with  the  interior  of  the  n-gon,  cuts  it  into  two  bounded  convex 
polygons.  The  number  of  sides  of  the  two  polygons  mi  and  m2  satisfy  the  following 
relation:  m\  + m2  = n + 4. 

Proof:  The  convexity  of  the  two  polygons  is  guaranteed  by  the  argument  that 
they  are  the  intersection  of  two  convex  sets  (the  original  convex  n-gon  and  the  two 
semi-planes,  from  Property  5.1  a semi-plane  is  a convex  set).  The  boundedness 
is  guaranteed  by  the  fact  that  they  both  are  subsets  of  the  original  n— gon.  It  only 
remains  to  show  that  the  relation  between  the  number  of  the  sides  of  the  two  polygons 
which  was  given  in  the  statement  of  Lemma  5.1  is  indeed  correct.  Suppose  that  one 
semi-plane  of  line  Ln  contains  k vertices  of  the  n-gon  and  the  other  semi-plane 
contains  n — k vertices,  where  1 < k < n — 1.  Obviously,  the  n-gon  is  cut  into  two 
polygons,  one  (&  + 2)-gon  and  one  (n  — k + 2)-gon  (see  Figure  5.2).  So  the  sum  of  the 
number  of  vertices  of  the  two  polygons  is:  mi  + m2  = (k  + 2)  + (n  — k + 2)  = n + 4. 
□ 

Lemma  5.2  Given  an  unbounded  convex  n-gon,  a straight  line  Ln  which  has  a com- 
mon point  with  its  interior  either: 

• cuts  it  into  an  unbounded  convex  m\-gon  and  a bounded  convex  m^-gon,  that 
satisfy  the  relation:  m\  + m2  = n + 4,  or 

• cuts  it  into  two  unbounded  convex  polygons  which  satisfy  the  relation:  mi+m2  = 
n + 3. 

Proof:  The  two  cases  are  shown  in  Figure  5.3  (lines  Ln(i ) and  Ln(ii)).  The  proof  is 
in  the  same  spirit  as  in  Lemma  5.1.  □ 
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Lemma  5.3  Suppose  a line  Ln  cuts  a convex  polygon  into  two.  The  two  new  polygons 
that  are  created  have  m\  and  m2  sides.  The  number  of  sides  of  the  new  polygons 
satisfy  the  following  relation:  max(mi,m2)  = n + 1. 

Proof:  Obvious.  If  one  notices  that  the  smallest  possible  number  of  vertices  for  a 
bounded  convex  n-gon  is  three  and  the  unbounded  convex  n-gon  with  the  smallest 
number  of  vertices  is  the  2-gon.  The  results  of  Lemmas  5.1  and  5.2  can  be  used  to 
show  that:  max(mi,m2)  = n + 1.  □ 

Remark:  In  Lemmas  5.1,  5.2  and  5.3  we  did  not  address  the  special  case  of  a 
line  passing  thourgh  a vertex  of  the  polygon.  Extension  of  our  results  to  these  cases 
is  straightforward.  In  the  proof  of  the  theorem  that  follows  we  will  not  consider  the 
case  when  a straight  line  is  parallel  with  one  of  the  edges  of  the  polygon.  Extension 
of  our  results  for  this  special  case  is  also  straightforward. 

Theorem  5.3  If  a set  of  n > 5 lines  defines  a convex  n-gon  this  n-gon  is  unique.  If 
a set  of  n > 3 lines  defines  a bounded  convex  n-gon  this  bounded  n-gon  is  unique. 

Proof:  Three  straight  lines  define  only  one  bounded  convex  3-gon  (triangle). 
Let  f0(x,y)  = 0,  f\(x,y)  = 0 and  f2(x,y)  = 0 be  the  equations  for  lines  L0,  L\  and 
Z/2  and  A,  B , C the  vertices  of  the  triangle  (see  Figure  5.4).  All  points  of  the  triangle 
satisfy  the  following  equations  [46,  page  134]: 

fi(x,y)  >0  for  i = 0,1,2.  (5.1) 

If  this  is  not  the  case,  one  can  always  substitute  fi(x,y ) = —fi(x,y)  to  make  it 
happen.  These  three  lines  also  define  three  unbounded  3-gons  which  satisfy  the 
equations  [46,  page  134]: 


fi(x,y)  > 0,  fi+i(x,  y)  > 0,  and  fi+2(x,y)  < 0 for  i = 0,1,2. 


(5.2) 
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Whenever  the  subscript  k , indicating  a member  of  a set  of  n elements,  exceeds  the 
largest  value  an  element  can  take  (i.e.,  n — 1 for  a set  of  n elements  when  counting 
starts  from  zero),  the  element  with  subscript  kmod{n ) will  be  implied.  Suppose  now, 
that  a fourth  line  L3  with  equation  f3(x,y)  = 0 is  added  to  the  set.  Two  cases  exist: 
either  L3  cuts  the  triangle  or  it  does  not  cut  it  (see  Figure  5.4).  It  can  be  proved  that 
the  two  cases  are  equivalent  by  just  reordering  the  set  of  lines  (if  the  line  L3  does  not 
cut  the  triangle  defined  by  lines  L0LiL2  then  the  line  L2  cuts  the  triangle  defined  by 
the  lines  L0L\L3,  this  fact  is  illustrated  in  5.4).  Hence,  we  only  need  to  check  what 
happens  to  the  unbounded  3-gons  when  a straight  line  cuts  the  triangle. 

Suppose  that  > 0 and  /3(C)  < 0 (see  Figure  5.5).  Obviously  L3  has 

no  common  points  with  the  boundary  of  the  unbounded  3-gon  L2BAL0  (it  intersects 
L0  and  L2  at  points  where  /1  > 0 and  E and  D belong  to  the  triangle  ABC)  and 
since  f3(A,B)  > 0 it  does  not  intersect  the  segment  AB.  Also,  for  the  unbounded 
3-gons  LqC BL\  and  L\ACL2,  it  is  obvious  that  they  have  at  least  one  common  point 
with  line  L3.  Because  line  L3  is  not  parallel  with  L\,  it  has  a common  point  with 
L\,  either  in  the  ray  AL\  or  in  the  ray  BL\.  So  a line  L3  that  cuts  the  triangle 
defined  by  lines  LqL\L2  always  gives  an  unbounded  convex  4-gon  by  cutting  one  of 
the  unbounded  3-gons  created  by  the  same  lines;  has  one  common  point  with  the 
second  unbounded  3-gon  giving  two  unbounded  3-gons;  and  has  no  common  point 
with  the  third  unbounded  3-gon.  This  means  any  set  of  four  straight  lines  gives  one 
bounded  convex  4-gon  and  one  unbounded  convex  4-gon  as  in  Figure  5.5. 

We  will  show  now  that  five  straight  lines  can  only  give  one  convex  5-gon,  either 
bounded  or  unbounded.  The  bounded  4— gon  is  determined  by  the  inequalities  /,•  > 0 
for  i = 0,1,2, 3 and  the  unbounded  4-gon  by  the  inequalities  /,  > 0 for  i = 0, 1 and 
fi  < 0 for  i = 2, 3.  If  this  is  not  the  case  we  can  always  make  it  happen  by  rearranging 
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(a) 


(b) 

Figure  5.4.  Cases  shown  in  Figures  (a)  and  (b)  are  equivalent.  If  a straight  line  Z3 
does  not  cut  the  triangle  formed  by  straight  lines  L0,  L\  and  L2 , one  of  the  lines  (L0, 
L\  or  L2  will  cut  the  triangle  formed  by  the  other  two  and  line  Z3. 
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Figure  5.5.  Line  L3  cuts  the  triangle  ABC  giving  a convex  bounded  4-gon  (ABED) 
while  cutting  the  unbounded  3-gon  LqCBL4  giving  an  unbounded  convex  4-gon 
(. LoCEFh ). 

the  set  of  four  lines.  If  line  L4  with  equation  f4(x,y)  — 0 cuts  the  4— gon  giving  a 
5-gon  and  a triangle,  one  of  the  following  has  to  be  satisfied: 

(i)  f4 (B)  < 0 ,f4(E,D,A)  > 0 (see  Figure  5.6,  Line  L4(i)), 

(ii)  /4(A)  < 0,  f4(B,  E,  D)  > 0 (see  Figure  5.6,  Line  L4(ii)), 

(in)  f*(D)  < 0,  f4(A,  B,  E)  > 0 (see  Figure  5.6,  Line  L4(iii)), 

(iv)  fi{E)  < 0,  f4(D,  A,  B)  > 0 (see  Figure  5.6,  Line  L4(iv)). 

We  will  show  that  if  one  of  these  four  cases  arises  there  is  no  way  the  unbounded 
4-gon  will  be  cut  by  the  straight  line  L4  to  give  a 5-gon.  First  we  study  case  (i) 
which  is  equivalent  with  case  (iii).  Because  f4(B)  < 0 and  f4(A)  > 0,  we  have 
f4(F ) < 0.  Hence,  line  L4  crosses  line  L3  within  the  segment  FE  (see  Lemmas  5.2 
and  5.3).  The  only  way  the  crossing  of  line  L4  to  the  unbounded  4— gon  LqCEFL4 
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Figure  5.6.  Cases  (i)-(iv)  illustrate  the  necessary  condition  for  a set  of  5 straight  lines 
to  define  a bounded  convex  5-gon. 

can  produce  a 5— gon  is  by  making  a triangle  one  of  whose  vertices  is  either  E or  F. 
This  implies  that  L4  has  a common  point  with  either  ray  FL\  or  segment  EC.  Both 
cases  are  impossible  because  L4  already  has  a common  point  with  the  line  L\ (within 
the  segment  AB)  and  line  L2  (within  the  segment  BE). 

Now  we  examine  case  (iv).  Because  f4(E)  < 0 and  f4(D ) > 0 we  have  f4(F)  < 0. 
Using  the  same  reasoning  we  can  show  that  /4(C)  < 0.  So  line  L4  crosses  lines  Lq 
and  L\  within  the  segments  CD  and  f?F,  respectively.  Also  line  L4  crosses  lines 
L2  and  L3  within  the  segments  EB  and  ED , respectively.  Obviously,  none  of  the 
points  of  intersection  of  L4  with  lines  T,,  i — 0, 1,2,3  belongs  to  the  boundary  of 
the  unbounded  convex  4— gon  LqC EF L\.  So  line  L4  has  no  common  points  with  the 
unbounded  convex  4-gon.  Case  (ii)  can  be  studied  in  the  same  way  as  case  (iv). 

Hence,  if  a line  L4  cuts  a bounded  convex  4-gon,  defined  by  a set  of  four  lines 
giving  a bounded  convex  5-gon,  it  is  guaranteed  that  the  unbounded  4-gon  which  is 
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defined  by  the  same  set  of  lines  will  not  give  a 5-gon.  We  can  also  prove  the  reverse, 
if  the  line  cuts  the  unbounded  convex  4-gon  defined  by  the  same  set  of  four  lines 
giving  a 5-gon,  it  will  not  cut  the  bounded  4-gon  defined  by  the  same  set  of  lines  to 
give  a 5-gon. 

We  proved  that  a set  of  5 lines  can  only  define  one  convex  5-gon,  either  bounded 
or  unbounded.  With  the  use  of  induction  we  will  prove  that  for  n > 5 a set  of  n lines 
can  define  only  one  convex  n— gon. 

1.  the  case  n = 5 has  already  been  proved, 

2.  suppose  that  if  a set  of  n = k > 5 lines  defines  a convex  k-gon,  this  k-gon  is 
unique, 

3.  we  will  prove  that  this  is  also  the  case  for  n = k + 1.  In  other  words  we  need 

to  prove  that  given  that  a set  of  lines  L0 , L\, . . . , Lk  defines  a (k  + l)-gon  this 
(k  + l)-gon  is  unique.  Obviously  any  subset  of  k lines,  say  L0,  L\, . . . , Lk_i, 
defines  at  least  one  convex  k— gon  (bounded  or  unbounded).  If  this  is  not  true 
then  because  of  Lemma  5.3  it  would  be  impossible  for  the  set  of  (k  + 1)  lines  to 
define  an  (k  + l)-gon.  By  (2)  the  k- gon  is  unique.  The  addition  of  one  more 
line  cuts  it  into  two  polygons  one  of  which  is  by  hypothesis  a (k+  l)-gon.  This 
( k + l)-gon  is  unique  because  the  k-gon  is  unique.  □ 

We  can  summarize  the  results  of  Theorem  5.3  as  follows  (the  assumption  is  that 
no  three  of  the  straight  lines  are  parallel  or  have  a common  point): 

1.  three  straight  lines  define  one  bounded  convex  3-gon  and  three  unbounded 


convex  3-gons, 
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Figure  5.7.  Figure  showing  the  case  where  two  lines  are  parallel,  (a)  The  four  lines 
define  a bounded  convex  4-gon;  (b)  the  four  lines  define  two  unbounded  convex 
4-gons. 

2.  four  straight  lines  define  one  bounded  convex  4-gon  and  one  unbounded  convex 
4-gon,  (if  two  of  the  lines  are  parallel  they  define  either  two  unbounded  convex 
4-gons  or  one  bounded  convex  4-gon,  see  Figure  5.7), 

3.  n straight  lines,  where  n > 5,  can  only  define  one  convex  rc-gon  (bounded  or 
unbounded). 

The  previous  theorem  helps  us  to  establish  the  following  propositions. 

Proposition  5.1  Any  bounded  convex  polygon  is  uniquely  determined  by  the  peaks  of 
its  Hough  transform. 

Proof:  Obviously  the  n sides  of  the  n— gon  give  n peaks  in  the  Hough  space.  By 
hypothesis,  these  n peaks  define  a bounded  convex  n— gon.  From  Theorem  5.3  this 
n— gon  is  unique.  □ 


Proposition  5.2  For  n > 5,  a convex  n—gon  is  uniquely  defined  by  the  peaks  of  its 
Hough  transform. 
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Figure  5.8.  An  example  where  two  non-convex  polygons  are  defined  by  the  same 
set  of  straight  lines.  Both  non-convex  hexagons  ADF El J A and  KECDGJK  are 
defined  by  the  same  set  of  straight  lines  (and  therefore  give  the  same  peaks  in  the 
Hough  space). 

Proof:  The  proof  is  similar  to  the  proof  of  Proposition  5.1.  □ 

5.3.2  Discussion 


We  showed  that  any  bounded  convex  polygon  is  uniquely  determined  by  the 
peaks  its  vertices  give  in  the  Hough  space.  We  also  showed  that  a convex  polygon 
with  more  than  five  sides  is  uniquely  determined  by  the  peaks  it  gives  in  the  Hough 
space.  Unfortunately,  as  one  can  conclude  by  the  example  of  Figure  5.8,  this  is  not 
the  case  for  non-convex  polygons.  In  the  section  that  follows  we  will  propose  an 
efficient  scheme  that  constructs  a convex  polygon  by  the  equation  of  its  edges. 
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5.4  An  Algorithm  for  the  Reconstruction  of  Convex  Polygons 

5.4.1  Motivation 

A problem  that  arose  in  the  previous  chapter  is  the  following:  if  we  have  seg- 
mented a set  of  n straight  lines  can  we  devise  an  efficient  algorithm  that  will  answer 
the  question  if  these  lines  define  a convex  n-gon?  If  the  answer  is  yes,  can  the  same 
algorithm  construct  the  convex  n-gon? 

5.4.2  Theory 

A precise  statement  of  the  problem  we  are  dealing  with  is  as  follows: 

Problem  5.1  Given  a set  of  n straight  lines  Lfii  = 0, . . . , n — 1)  find  if  they  define  a 
convex  n-gon  and,  if  it  exists,  construct  it. 

The  algorithms  presented  in  this  section  make  use  of  the  following  Lemmas. 

Lemma  5-4  A necessary  and  sufficient  condition  for  a straight  line  Ln  to  cut  a bounded 
convex  n-gon  giving  an  (n  + \)-gon  is  that  for  one  vertex  of  the  n—gon: 

fn(Pi)  < 0 and  fn(Pi_1,Pi+ 1)  > 0.  (5.3) 

Proof:  Necessary:  Obvious.  Since  the  straight  line  Ln  cuts  the  polygon  giving  an  (n  + 
1)— gon,  one  of  the  semi-planes  it  defines  contains  n — 1 vertices  of  the  polygon  and  the 
other  contains  1 vertex.  Therefore,  fn(Pi)  < 0 and  fn(P0,  • • • , Pi-u  Pi+i,  • • • , Pn-\)  > 
0 which  implies  fn(Pl)  < 0 and  fn(Pt_uPi+1)  < 0. 

Sufficient:  Line  Ln  cannot  have  more  than  two  common  points  with  the  boundary  of 
the  convex  set  (Theorem  5.2).  It  already  has  one  common  point  with  the  edge  Pi-\Pi 
and  one  with  the  edge  P,Pt+i . Therefore  it  does  not  intersect  any  other  edge  of  the 
bounded  convex  polygon.  This  means  that  edges  P0  . . . Pi-\Pi+\  . . . P„_i  belong  to 
the  same  semi-plane  and  therefore  a bounded  convex  (n  + l)-gon  is  defined  by  the 
intersection  of  that  semi-plane  and  the  bounded  convex  polygon.  □ 
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Lemma  5.5  Given  an  unbounded  convex  n-gon  LqPq  . . . Pn-2Ln-\  and  a straight  line 
Ln  with  equation  fn(x,y)  = 0.  Given  that  V\  is  the  point  of  intersection  of  lines  Lq 
and  Ln  and  V2  o/Tn_i  and  Ln,  a necessary  and  sufficient  condition  for  straight  line 
Ln  to  cut  the  n-gon  giving  an  ( n + l)-gon  is  that  one  of  the  following  is  satisfied: 

1.  for  a vertex  Pi  with  i = 1, . . . , n — 3,  fn{Pi ) < 0 and  /n(P;_i,  Pi+i)  > 0, 

2.  fn(Po)  < 0,  fn{Pi)  > 0 and  fi{\\)  > 0, 

3.  fn(Pn-2 ) < 0,  fn{Pn-z)  > 0 and  fn-2(V i)  > 0, 

l fn(Po)  > 0,/n(Pn_2)  > 0 and  one  of  the  following 

(a)  fi(Vi)  > 0, 

(b)  fn-2(V2)  > 0, 

(c)  fi{Vi)  > 0 and  fn-2(V 2)  > 0. 

Proof:  Omitted.  □ 

Lemmas  5.4  and  5.5  can  be  used  to  propose  an  0(n2)  algorithm  for  the  solution 
of  the  problem.  The  checking  routine  of  such  an  algorithm  for  the  bounded  convex 
polygon  case  would  be: 
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Algorithm  f 0(n2)  for  construction  of  a convex  polygon  by  its  edges 

begin:  check  (bounded  convex  n— gon  (PqP\  ■ ■ ■ Pn- 1),  line  Ln ) 
for  i = 0 to  n — 1 

if  fn(Pi)fn(Pi-i)  < 0 and  fn(Pt)fn(Pl+ 1)  < 0 
find  V\  G LnC\  Pi  Pi- 1 and  Vi  € Ln  D PiPi+i 

return(bounded  convex  (n  + l)-gon  P0Pi  . . . Pi-i Vi V^-Pi+i  . . . Pn-\) 

end 

end 

return  (“NO  BOUNDED  CONVEX  POLYGON  DEFINED”) 
end:check 

end 

What  this  algorithm  does  is  check  if  the  vertices  of  the  convex  polygon  satisfy 
the  conditions  of  Lemma  5.4.  If  this  is  the  case  for  a vertex  P,-,  it  finds  the  two 
new  vertices  Vi  and  V2  and  replaces  Pi  with  V1V2.  After  doing  the  replacement  it 
returns  the  new  polygon.  If  no  substitution  takes  place  it  continues  to  the  next  vertex 
until  no  vertex  remains.  If  no  vertex  satisfies  the  conditions  of  Lemma  5.4,  a report  is 
issued  that  this  set  of  n lines  does  not  define  a convex  n-gon.  A similar,  slightly  more 
complicated,  routine  can  be  used  for  unbounded  convex  polygons.  If  one  builds  the 
initial  bounded  and  unbounded  3-gon  and  proceeds  incrementally  (one  straight  line 
at  a time)  the  algorithm  would  require  0(n2)  computations.  We  state  two  theorems 
that  will  help  us  propose  an  O(nlogn)  algorithm. 
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Figure  5.9.  A bounded  convex  6-gon.  It  is  easy  to  see  that  C2  > C\  > Co  and 
c5  > c4  > c3  and  Co  < C5  and  C3  < C2. 

Theorem  5.1  A bounded  convex  n-gon  can  always  be  cut  into  two  pieces  LqL\  ...  Li 
and  Li+xLi+2  . . . A„_i  (as  we  travel  it  counterclockwise)  such  that  for  each  piece  ci  < 
C/_|_2  j and  Ci  ^ c,_ 1-1 ; cn — 1 cq. 

« 

Proof:  Omitted.  □ 

Figure  5.9  can  be  used  to  visualize  what  is  stated  in  Theorem  5.4.  The  polygon 
of  Figure  5.9  can  be  cut  into  two  pieces,  namely:  L0L1L2  and  L3L4L5.  We  also  need 
to  notice  that  Theorem  5.4  is  equivalent  to  Theorem  3.6  of  Preparata  and  Shamos 
[69,  page  99].  We  state  the  equivalent  of  Theorem  5.4  for  unbounded  convex  (UC) 
n-gons. 

Theorem  5.5  For  an  unbounded  convex  n-gon  LqL\  . . . T„_i  one  of  the  following  is 
satisfied: 


1.  ci  < c;+ 1 for  0 < l < n — 1, 
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2.  it  can  be  cut  into  two  pieces,  LqL\  . . . Li  and  Li+\Li+2  • • ■ Ln-\  such  that  for 
each  piece  ci  < q+i  and  c;  > Cj+i . 

Proof:  Omitted.  □ 

Figures  5.10(a)  and  5.10(b)  illustrate  what  is  stated  in  Theorem  5.5.  The  un- 
bounded convex  polygon  of  Figure  5.10(a)  satisfies  Condition  1 of  Theorem  5.5  while 
the  unbounded  convex  6-gon  of  Figure  5.10(b)  satisfies  Condition  2 of  Theorem  5.5 
(it  can  be  cut  into  two  pieces  L0L\ . . . L4  and  L5).  The  following  definition  will  be 
useful  in  the  theory  that  will  follow: 

Definition  5.1  The  vertices  of  a convex  polygon  L0L\  . . . Ln-\,  Pi  (points  of  intersub- 
section of  line  Li-i  and  Li)  where  c,_i  > c,-  will  be  called  Orientation  Change  Vertices 
(OCV). 

If  i > n — 1 we  take  imod(n).  Obviously  a bounded  convex  (BC)  n— gon  always 
has  two  OCV  and  an  unbounded  convex  (UC)  n-gon  can  not  have  more  than  one 
OCV.  Lemma  5.4  and  Theorem  5.4  can  be  combined  to  give  the  following  Lemma: 

Lemma  5.6  Given  the  BC  n-gon  defined  by  the  lines  LqLi  . . . T„_i  and  a straight 
line  Ln  with  cn  > max(co , . . . , cn_ i);  a necessary  and  sufficient  condition  for  the  line 
Ln  to  cut  the  n-gon  giving  a BC  (n  + l)-gon  is  that  if  the  Pq  and  Pi  are  the  OCVs 
of  the  n—gon  either: 

L fn(Po)fn(Pn-i)  < 0 and  fn{Po)fn(Pi)  < 0,  or 

fn{Pi)fn(Pi- 1)  < 0 and  fn(Pi)fn(Pi+ 1)  < 0. 

Proof:  Omitted.  ED 

Figure  5.11  shows  the  two  cases  of  Lemma  5.6  for  the  bounded  convex  polygon 
of  Figure  5.9. 
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Figure  5.10.  (a)  An  unbounded  convex  5-gon  that  has  no  Orientation  Change  Vertex, 
(b)  An  unbounded  convex  6-gon  with  an  Orientation  Change  Vertex  ( P4 ). 
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Figure  5.11.  Line  Lq{i)  satisfies  the  conditions  of  Lemma  5.6  (Case  1)  and  Line  L6(n) 
satisfies  the  conditions  of  Lemma  5.6  (Case  2 with  i = 3). 
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In  a similar  fashion  Lemma  5.5  and  Theorem  5.5  can  be  combined  to  give  a 
condition  for  the  unbounded  convex  n- gon  case. 

Lemma  5.7  Given  an  unbounded  convex  n-gon  defined  by  the  lines  LoL\ . . . Ln_\  and 
a straight  line  Ln  with  cn  > max(c0, . . . , c„_i).  A necessary  and  sufficient  condition 
for  the  line  Ln  to  cut  the  n—gon  giving  a convex  (n  + 1 )-gon  is  that  either: 

1.  the  n-gon  does  not  have  an  OCV,  and  V\  € Ln  fl  Lq  and  V2  € Ln  fl  Ln_i  and 
either: 

(a)  fi(V\)  > 0 giving  an  UC  (n  + l)-gon  with  an  OCV  (V\),  or 

(b)  y”n_2  ( V2 ) > 0 giving  an  UC  (n  + 1 )-gon  without  an  OCV. 

2.  the  n—gon  has  an  OCV  Pi,  and  V\  € Ln  fl  Li,  and  and  V2  € Ln  fl  L,+ 1,  then 

(a)  if  i > 0,  and  i < n — 2 then  fn(Pi ) < 0 and  fn(Pi- 1,  P+i)  > 0, 

(b)  ifi  = 0 then  fn(P0 ) < 0,  fn{P\)  > 0 and  fi(Vi)  > 0 

(c)  ifi  = n - 2 then  fn(Pn-2 ) < 0,  /„(T,„-3)  > 0 and  fn-2(V 2)  > 0 

(d)  if  fn(Po,  Pi  1 Pn-2)  > 0 then  fi{Vi)  > 0 and  /„-2(Vr2)  > 0 

Proof:  Omitted.  □ 

Figure  5.12  illustrates  the  cases  of  Lemma  5.7. 

5.4.3  The  algorithm 

Below  we  provide  some  details  for  each  part  of  the  algorithm  that,  given  a set 
of  n straight  lines,  constructs  the  n-gon  they  define,  after  checking  if  it  exists. 

Part  1 

In  this  part  sorting  of  the  straight  lines  with  increasing  slope  takes  place.  Sorting 
algorithms  like  the  heapsort  or  the  mergesort  [14,  76]  can  be  used.  It  is  well  known 
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(a) 


(b) 


(c) 

Figure  5.12.  Figure  illustrating  the  conditions  of  Lemma  5.7.  (a)  Lines  Ls(f,  ii ) satisfy 
the  conditions  of  Lemma  5.7  (Cases  1(a)  and  1(b)  respectively);  (b)  line  £7(2)  satisfies 
the  conditions  of  Lemma  5.7  (Case  2(a));  (c)  lines  L6(i,ii ) satisfy  the  conditions  of 
Lemma  5.7  (Cases  2(b)  and  2(c)). 
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that  given  a set  of  n straight  lines  the  computation  time  for  these  algorithms  will  be 
O[nlogn). 

Part  2 

This  part  takes  as  input  the  three  lines  with  the  smallest  c’s  and  gives  as  output 
the  four  convex  polygons  (one  triangle  and  three  unbounded  3-gons)  defined  by  these 
three  lines.  Given  that  cq  < c\  < c2  it  is  easy  to  show  that  the  output  will  always  be 
(see  Figure  5.13): 

1.  bounded  convex  3-gon  (triangle)  LXL0L2  (with  OCV’s  L\  fl  L0  and  Lx  fl  T2), 

2.  unbounded  convex  3-gon  L2L0LX  (with  OCV  L2  fl  T0), 

3.  unbounded  convex  3-gon  LxL2Lq  (with  OCV  L2  fl  To),  and 

4.  unbounded  convex  3-gon  L0LXL2  (without  an  OCV). 

This  part  of  the  algorithm  takes  0(1)  time. 

Part  3 

In  this  part  Lemmas  5.6  and  5.7  are  used  iteratively  to  construct  the  new  poly- 
gon, after  we  check  if  it  exists.  Since  we  need  to  consider  only  5-6  cases  per  iteration 
this  part  of  the  algorithm  needs  0(n ) computations. 

5.4.4  Examples 

In  this  Section  we  will  present  two  examples  that  illustrate  how  the  algorithm 


works. 
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Figure  5.13.  Three  straight  lines  defining  four  convex  3-gons  (one  bounded  and 
three  unbounded).  Notice  that  if  Co  < c\  < c2  the  counterclockwise  order  for  lines  of 
the  bounded  3-gon  is:  L\L0L2  while  for  the  unbounded  convex  3-gons  the  order  is: 
L2LqL\,  L\L2Lq  and  LqL\L2. 


Example  1 


A set  of  seven  straight  lines  is  given.  The  parameters  of  these  lines  are: 


L = {(c,-,ra;),  for  i — 0 to  6} 


-11.5087 

0.8420 

-5.5757 

5.6243 

1.0975 

-0.6976 

-25.4152 

24.0836 

-0.0412 

0.0754 

-0.4391 

0.7500 

0.0282 

0.9676 

(5.4) 


After  the  first  step  (sorting  with  respect  to  c)  the  parameters  that  are  input  to 


-25.4152 

24.0836 

-11.5087 

0.8420 

-5.5757 

5.6243 

-0.4391 

0.7500 

-0.0412 

0.0754 

0.0282 

0.9676 

1.0975 

-0.6976 

the  second  step  are: 
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The  output  of  the  second  part  consists  of  four  3— gons  (see  Figure  5.14(a)): 

1.  triangle  ABC, 

2.  unbounded  3-gon  L2BALi  with  OCV  B, 

3.  unbounded  3-gon  L\CBLo  with  OCV  B , and 

4.  unbounded  3-gon  LoACL2  without  an  OCV, 

where  A G L0  0 L\,  B £ L0  0 L2,  and  C € L\  0 L2.  The  coordinates  of  these  points 
are:  A(1.6713, -18.3923),  £(0.9304,0.4364)  and  C(-0.8061, 10.1187).  In  the  third 
part  of  the  algorithm  we  follow  the  four  convex  polygons  (i)-(iv)  individually. 

First  we  check  the  triangle  ABC.  The  OCVs  for  this  bounded  convex  polygon 
are  the  points  A and  C.  For  the  line  L 3 we  get  f 3(A) f 3(B)  < 0 and  /3(A)  f3(C)  < 0. 
Since  E e Lx  C i3  and  D € L0  0 L3  with  £(0.0668, 0.0727)  and  £(0.9462,  0.0364),  D 
substitutes  for  A as  OCV  (see  Figure  5.14(b)). 

For  the  line  L4  we  get  f4(C)f4(E)  < 0 and  /4(C)/4(£)  < 0.  We  calculate  the 
coordinates  of  points  F 6 L2  fl  L4  and  G £ L\  fl  L4.  We  find  that  £(0.8310,0.9910) 
and  G(— 0.0109, 0.9673).  Point  G substitutes  for  C as  OCV  (see  Figure  5.14(c)). 

The  pentagon’s  vertices  are  (in  a counterclockwise  order):  DBFGE.  The  OCVs 
for  this  polygon  are  vertices  D and  G.  Therefore,  when  line  L5  is  added  we  only  need 
to  check  what  happens  to  these  two  vertices.  Proceeding  in  the  same  way,  we  see 
that  /5(G)/5(£)  < 0 and  fs(G)fs(E)  < 0 and  / is  the  new  OCV  of  the  polygon 
(see  Figure  5.14(d)).  The  coordinates  of  vertices  H and  / are:  £(0.5297,0.9826) 
and  7(0.0072,0.7534).  Repeating  this  process  for  line  L&  we  get  the  bounded  con- 
vex polygon  KBFHIEJ  (see  Figure  5.14(e)).  The  coordinates  of  J and  K are: 
J(0. 6789, 0.0475)  and  I\  (0.9347, 0.3282).  We  would  have  proceeded  in  the  same  man- 
ner if  more  than  seven  lines  belonged  to  the  set  that  was  input  to  our  algorithm.  In 
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the  same  way  we  check  the  unbounded  convex  3-gon  defined  by  the  first  three  straight 
lines  and  we  find  that  no  other  7-gon  is  defined  by  these  lines. 

Example  2 


In  this  example  we  are  dealing  with  an  unbounded  convex  n-gon.  Unbounded 
convex  polygons  are  described  by  a set  of  pairs  of  coordinates.  The  first  and  the  last 
pair  of  this  set  are  the  parameters  of  the  bounding  rays.  This  is  consistent  with  the 
notation  discussed  in  Section  5.2. 

The  input  is  a set  of  six  straight  lines  with  the  following  parameters: 


L = {(c;,ra8),  for  i — 0 to  5} 


' -6.4430  1.0111 

-1.3120  0.5958 

—0.1126  0.3326 

0.6635  -0.0599 

4.7861  -3.2019 

15.2375  -11.7947 


(5.5) 


After  the  second  part  we  find  that  the  first  three 
unbounded  convex  3— gon  (see  Figure  5.15(a)): 


lines  define  the  following 


‘ -6.443  1.011  ‘ 

0.08083  0.4897 
0.2194  0.3079 

-0.1126  0.3326 


(5.6) 


This  is  an  unbounded  convex  polygon  without  an  OCV  (the  first  and  last  members  of 
the  set  CP  are  the  bounding  rays  of  the  unbounded  convex  polygon,  while  the  second 
and  third  members  are  the  vertices  of  that  polygon  in  counterclockwise  order).  The 
fourth  straight  line  of  the  set  satisfies  the  conditions  of  Lemma  5.7  (Case  1(b))  with 
Vi  = C = (0.5058,0.2757)  (see  Figure  5.15(b)).  Hence,  the  new  unbounded  convex 
polygon  is: 


-6.443 

1.011 

0.08083 

0.4897 

0.2194 

0.3079 

0.5058 

0.2757 

0.6635 

-0.0599 

CP  = 


(5.7) 
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(<0  (d) 


(e) 


Figure  5.14.  Figure  clarifying  the  steps  of  Example  1. 


124 


Proceeding  the  same  way  we  can  find  that  the  fifth  line  satisfies  Lemma  5.7  (Case 
1(b)),  with  Vi  = D — (0.7621,0.4457)  (see  Figure  5.15(c)).  The  new  unbounded 
convex  polygon  is: 


CP 


is: 


CP 


-6.443 

1.011 

0.08083 

0.4897 

0.2194 

0.3079 

0.5058 

0.2757 

0.7621 

0.4457 

4.7861 

-3.2019 

ee  Figure 

5.15(d))  tl 

-6.443 

1.011 

0.08083 

0.4897 

0.2194 

0.3079 

0.5058 

0.2757 

0.7621 

0.4457 

0.8222 

0.7331 

15.2375 

-11.7947 

(5.8) 


(5.9) 


5.4.5  Results-Conclusions 


The  procedure  for  verification  of  the  correctness  of  the  algorithm  is  as  follows:  a 
random  sample  of  n points  within  the  region  (0, 1)  x (0, 1)  was  chosen.  Using  available 
software  the  convex  hull  of  these  points  was  found.  After  that,  the  equations  of  the 
edges  were  derived  and  rearranged  randomly.  The  set  of  equations  of  these  straight 
lines  was  the  input  to  our  code.  The  differences  of  the  coordinates  of  the  vertices 
of  the  convex  hull  as  provided  by  the  convex  hull  building  program  and  our  code 
were  compared  and,  except  for  round  off  errors,  no  differences  were  found  for  thirty 


thousand  cases. 
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(c) 


(d) 


Figure  5.15.  Figure  clarifying  the  steps  of  Example  2. 


CHAPTER  6 

CONCLUSIONS-FUTURE  WORK 


In  this  chapter  we  review  the  results  of  our  work  and  propose  future  research 
directions. 

6.1  Contributions  of  this  Thesis 
6.1.1  Hough  Transform  Understanding 

We  used  Digital  Topology  theory  to  develop  an  understanding  of  the  factors 
that  influence  the  Hough  transform.  The  study  of  the  spreading  of  a DSL  segment, 
peak  extension,  choice  of  A p,  and  performance  prediction  provides  us  with  a solid 
understanding  of  the  effects  that  influence  the  accuracy  of  the  Hough  transform. 

Spreading 

A new  equation  that  predicts  the  worst  case  spreading  as  a function  of  the 
length  of  a DSL  segment  was  proved.  The  new  equation  differs  from  previous  attemps 
[88,  53]  because  it  treats  the  DSL  segments  as  digital  and  not  Euclidean  objects.  The 
use  of  Digital  Topology  theory  towards  this  goal  was  instrumental.  The  correctness 
of  the  formula  was  shown  experimentally. 

Peak  Extension 

It  was  proven  that  even  in  the  case  of  a DSL  segment  being  exactly  collinear 
with  a strip,  the  votes  it  gives  to  the  particular  strip  may  be  equal  to  the  votes  it 
gives  to  its  neighboring  strips  along  the  0-direction.  This  happens  when  the  DSL 
segment  totally  belongs  to  the  intersection  of  these  strips.  It  was  also  shown  that  the 
interplay  of  the  peak  extension  and  the  peak  spreading  may  lead  to  situtations  where 
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the  maximum  number  of  votes  is  attributed  to  a cell  that  does  not  best  represent  the 
parameters  of  the  DSL  segment. 

Choice  of  A p 

We  showed  that  if  we  want  the  members  of  the  parameter  space  to  represent 
DSLs  there  has  to  be  an  adaptive  A p = h max( | cos  6 | , |sin0|).  If  A p is  not  equal  to 
this  optimum  value  most  of  the  members  of  the  parameter  space  do  not  represent 
DSLs. 

Performance  Prediction 

The  use  of  Digital  Topology  theory  made  the  study  of  performance  prediction 
an  easy  task.  In  fact,  in  our  study  of  the  effects  of  spatially  uniform  noise,  we  were 
able  to  take  into  account  the  effects  of  the  digitization  of  the  image  and  the  parameter 
space.  Theoretical  predictions  generally  were  close  to  measured  values.  The  large 
discrepancies  that  appeared  in  some  directions  were  attributed  to  non-uniformities. 
Members  of  the  parameter  space  sometimes  contain  more  pixels  than  predicted  by 
theory.  This  is  the  result  of  a suboptimal  choice  of  A p = h.  If  A p is  chosen  according 
to  Digital  Topological  criteria  (see  Subsection  2.5)  these  effects  do  not  appear.  These 
effects  were  shown  to  become  important  in  cases  when  we  are  dealing  with  highly 
noisy  environments  or  environments  where  large  concentrations  of  feature  points  like 
texture  appear. 

Postprocessing 

A theoretical  framework  was  provided  for  the  calculation  of  the  length  of  DSL 
segments  from  the  Hough  transform.  The  proposed  algorithm  takes  into  account 
noise  and  other  objects  that  may  appear  in  the  image.  Local  uniformity  of  the  noise 
is  the  underlying  assumption  of  the  algorithm.  This  assumption  is  not  a very  strict 
one  and  it  almost  always  appears  in  practice.  We  showed  that  our  length  calculation 
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algorithm  is  equivalent  to  a filtering  of  the  Hough  space.  We  proved  that  the  use  of 
this  filtering  method  aggravates  the  noise,  but  if  we  assume  a uniform  background 
noise  we  will  get  a zero  mean  filtered  space.  Compared  to  other  filtering  techniques 
[55,  71]  our  filtering  scheme  has  the  distinct  advantage  of  an  intuitive  threshold 
selection. 

6.1.2  Object  Recognition/Representation 
Polygon  Recognition 

We  proposed  a Hough  transform-based  algorithm  for  the  recognition  of  convex 
polygons.  The  algorithm  is  data-driven  and,  therefore,  no  a priori  knowledge  of  the 
objects  to  be  recognized  is  needed.  The  algorithm  is  not  at  its  best  when  dealing 
with  complicated  environments  that  contain  many  polygonal  objects.  In  these  envi- 
ronments the  required  computation  increases  combinatorially  and  the  larger  number 
of  possible  matches  makes  more  likely  the  possibility  of  false  positives.  An  advantage 
of  the  method  is  the  capability  of  dealing  with  occlusions,  cluttered  environments 
and  the  computational  efficiency  when  compared  with  some  earlier  methods  (e.g., 
generalized  Hough  transform). 

Polygon  Representation/Reconstruction 

We  proved  that  a bounded  convex  polygon  is  uniquely  represented  by  the  param- 
eters of  its  vertices.  We  also  showed  that  this  is  the  case  for  all  convex  polygons  with 
more  than  four  sides.  We  provided  counter-examples  to  show  that  for  non-convex 
polygons  this  uniqueness  of  representation  does  not  hold.  Finally  we  proposed  an 
efficient  algorithm  for  the  reconstruction  of  a convex  polygon  from  the  parameters  of 
its  edges. 
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6.2  Future  Work 

Hough  transform  Understanding 

Even  though  we  improved  the  understanding  of  the  various  factors  and  their 
effects  on  the  Hough  transform  many  questions  still  remain  unanswered.  We  showed 
that  the  peak  cell  does  not  necessarily  best  represent  the  parameters  of  a DSL  seg- 
ment. This  is  the  case  whether  this  cell  is  in  the  peak  in  the  Hough  domain  or  the 
filtered  Hough  domain.  The  question  that  naturally  arises  is:  ‘is  there  any  way  we 
can  achieve  better  accuracy  in  the  estimated  parameters  from  the  Hough  transform 
space?’.  This  question,  essentially,  leads  to  the  proposal  of  another,  more  intelli- 
gent, criterion  for  “peak”  selection.  The  work  of  Niblack  and  Petkovic  [64]  is  an 
attempt  towards  this  goal  but  relies  on  the  unrealistic  assumption  of  the  absence  of 
noise.  Other  attemps  relied  on  multiresolution  representations  to  achieve  improved 
accuracy  [5,  72,  38]. 

Extensions  to  Radon  Transform 

Deans  has  shown  that  the  Hough  transform  is  a special  case  of  the  well-known 
Radon  transform  [20].  The  Radon  transform  has  found  applications  in  a number  of 
scientific  areas  such  as  astronomy  [73],  image  processing  [45],  and  medical  imaging 
[62,  67].  The  study  of  phenomena  similar  to  the  ones  we  studied  in  Chapter  2 for  the 
Radon  transform  would  be  very  interesting. 

Extensions  to  Other  Shapes 

A wealth  of  Hough  transform-based  shape  recognition  algorithms  has  appeared 
in  current  literature  [7,  19,  35,  47,  63].  The  study  of  the  effects  of  digitization  errors 
to  the  parameter  space  for  shapes  such  as  circles  and  ellipses  would  be  interesting. 
Hu  and  Destine  studied  the  effects  of  uniform  noise  on  the  parameter  space  for  circles 
[32].  They  assumed  that  both  image  and  parameter  spaces  are  continuous.  Efforts 
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should  be  directed  towards  the  extension  of  their  results  for  digitized  image  and 
parameter  spaces. 

Object  Recognition 

As  we  have  stated,  the  idea  of  incorporating  length  information  into  an  object 
recognition  scheme  has  been  around  for  some  time  [25].  On  the  other  hand  no  com- 
putational implementation  of  this  idea  has  been  presented.  All  model-based  methods 
utilize  the  number  of  votes  of  the  peaks.  We  have  shown,  that  due  to  the  spreading 
and  noise,  the  number  of  votes  of  the  members  of  the  parameter  space  does  not  pro- 
vide any  meaningful  information  about  the  length  of  the  DSL  segment.  Therefore, 
there  is  a need  to  incorporate  our  length  estimation  scheme  into  the  polyhedral  object 
recognition  algorithm. 

Extension  of  our  convex  polygon  recognition  scheme  to  generic  polygons  would 
also  be  important.  The  new  scheme  should  be  model-based  and  should  exploit  paral- 
lelism, collinearity  and  other  non-accidental  properties  [8,  58]  towards  efficient  group- 
ing of  the  primitives. 

Accuracy  of  the  Recognized  Objects 

An  important  issue  that  current  literature  has  not  dealt  with  is  the  quantitative 
assessement  of  the  accuracy  of  the  object  localization.  To  the  best  of  our  knowledge, 
no  Hough-based  recognition  scheme  has  provided  a theoretical  framework  for  the 
errors  of  the  localization  of  objects.  It  is  hoped  that  the  discussion  of  the  digitization 
errors  algong  with  the  calculations  of  Chapter  3 will  provide  the  motivation  for  a 
further  exploration  of  accuracy  issues  of  the  Hough  transform. 


APPENDIX  A 
PROOF  OF  EQUATION  3.4 


If  <p  is  the  angle  of  the  slope  of  the  DSL  segment,  then  for  (p  E [0, 7r/ 4) , from 
the  definition  of  the  (p,  #)  parameterization  of  the  line,  one  can  find  # = 7t/2  + <p. 
Supposing  that  the  vertices  of  a segment  have  coordinates:  (x\,y\)  and  (#2, 2/2)5  then 
cos  <p  = \x2  — xi\h/l  where  / is  the  length  of  the  DSL  segment,  |#2  — «i|  is  the  absolute 
value  of  x2  — X\  and  h the  size  of  the  pixel.  It  is  easy  to  see  that  for  (p  E [0, 7r/ 4) 
the  relation  \x2  — X\\  — vl  — \ holds  where  vl  is  the  number  of  pixels  of  the  DSL 
segment.  Therefore  the  equation  that  relates  the  number  of  pixels  that  belong  to  the 
segment,  with  its  length  is:  / = (vL~l)h/  cos  (p  — ( VL~l)h / cos(</>  — 7t/2).  Because  for 
6 E [7r/2,37r/4)  the  relation  cos(#  — 7t/2)  = sin#  holds,  we  have  l = ( vl  — 1 )h/  sin#. 

For  (j)  E [7t/4, 7r/2)  we  have  \y2  — yi\  = vl~  1.  Considering  that  for  # E [37t/4,7t), 
sin(#  — 7t/2)  = —cos#  = |cos#|,  we  can  show  that  / = ( vl  — l)h/s'm(f)  — ( vl  — 
1 )h/  sin(#  — 7t/2)  = ( vl  — l)/*/|cos#|.  If  we  carry  out  the  same  calculations  for  the 
rest  of  the  range  of  (f>  (two  more  cases,  (p  E [7t/2,  37t / 4)  and  (p  ^ [37r/4,  tt))  we  obtain 
the  following  generic  formula: 

{ = (vL~l)h 
max(|  cos  #|,  | sin  #|) 


(A.l) 
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APPENDIX  B 
STATISTICS  TESTS 


In  this  Appendix  we  present  our  experiments  that  confirm  our  claim  that  if  the 
noise  is  spatially  uniform  the  votes  of  a cell  (i,j)  of  the  parameter  space  follow  a 
binomial  distribution.  In  Section  B.l  we  present  the  theoretical  background  while  in 
Section  B.2  we  present  the  methodology  of  our  experiments. 

B.l  Goodness  of  fit  tests 

Goodness  of  fit  tests  [23]  are  used  to  determine  if  the  outcome  from  an  experi- 
ment fit  a certain  probability  distribution.  According  to  the  tests  the  null  hypothesis 
H0  (i.e.,  H0:  sample  B is  from  a distribution  A)  is  checked  against  an  alternative 
hypothesis  Ha  (i.e.,  Ha : sample  B does  not  come  from  distribution  A).  A significance 
level  a is  specified.  This  constant  quantifies  the  level  of  confidence  we  would  like  to 
have  on  the  outcome  of  our  data  processing.  The  sample  space  C is  partitioned  into 
k subintervals,  Co,  C\, . . . , Ck- 1 and  a test  statistic  y2  is  calculated  according  to  the 
following  formula: 

x2  = gfe^)!  (B.l) 

t=0  e* 

where  y,  is  the  measured  number  of  outcomes  for  the  category  C\,  et  = nP{C{), 
n = J2i=oVi-  The  null  hypothesis  is  accepted  if  y2  < Xa,w>  where  x2  „ is  the  value 
the  x2  distribution  function  takes  for  a level  of  confidence  and  v degrees  of  freedom. 
The  parameter  v is  calculated  from  the  following  equation: 

v = k — 1 — r (B.2) 

where  r is  the  number  of  parameters  of  the  distribution  that  were  calculated  from 
the  sample  and  k is  the  number  of  subintervals  used  for  the  calculation  of  x2- 
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B.2  Experiment 

For  a number  of  randomly  chosen  strips  we  proved  experimentally  that  their  dis- 
tribution is  a binomial  with  a 95%  level  of  confidence.  The  steps  of  the  experimental 
process  were  the  following: 

1.  Perform  the  Hough  Transform  for  an  image  whose  pixels  were  all  feature  points. 
With  this  method  we  were  able  to  identify  the  number  of  pixels  n belonging  to 
each  strip. 

2.  The  image  was  corrupted  with  5%  uniform  noise  and  the  number  of  votes  for 
the  strips  of  interest  was  reported.  This  process  was  repeated  for  a thousand 
times. 

3.  The  results  of  the  experiment  were  processed  according  to  the  Goodness  of  fit 
test  described  in  the  previous  section. 

(a)  Hq\  the  collected  votes  for  the  strip  come  from  a binomial  distribution, 
Ha:  the  collected  votes  for  the  strip  do  not  come  from  a binomial  distri- 
bution, 

(b)  level  of  confidence  chosen:  a = 0.95, 

(c)  in  partitioning  the  sample  space  we  followed  the  guidelines  provided  by 
Snedecor  and  Cochran  [83,  page  77].  We  can  calculate  the  test  statistic 
X2  from  Equation  B.l:  no  class  expectation  was  below  1 except  for  the 
two  extreme  expectations  and  most  of  the  other  expected  values  exceeded 
5.  The  values  of  ?/,•  and  are  plotted  in  Figures  B.l,  B.2,  and  B.3.  From 
Tables  [2]  we  found  X0.99  v and  compared  it  with  the  x2  we  computed. 

Experimental  results  for  three  members  of  the  parameter  space,  (i,j)  = {- 
(10, 22),  (55,  7),  (200, 22)},  are  reported  in  Tables  B.l,  B.2,  B.3.  The  number  of  points 
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Figure  B.l.  Comparison  of  theoretical  and  experimental  Probability  Density  Function 
(PDF)  for  the  case  (i,j)  = (10,22). 


Figure  B.2.  Comparison  of  theoretical  and  experimental  Probability  Density  Function 
(PDF)  for  the  case  (i,j)  = (55,7). 

belonging  to  these  strips  was  found  to  be:  154, 256  and  264  respectively.  These  Tables 
report  the  results  for  the  three  cases  shown  in  Figures  B.l  B.2,  and  B.3. 

In  Table  B.4  we  compare  the  y2  that  we  computed  for  the  three  cases  of  Tables 
B.l,  B.2  and  B.3.  The  degrees  of  freedom  v were  calculated  from  Equation  B.2  with 
r = 0 (we  used  our  a priori  knowledge  of  the  mean  value  instead  of  calculating  it 
from  the  sample).  Obviously  for  all  three  cases  the  null  hypothesis  is  accepted. 
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Figure  B.3.  Comparison  of  theoretical  and  experimental  Probability  Density  Function 
(PDF)  for  the  case  (i,j)  = (200,22). 


Table  B.l.  Goodness  of  fit  test  data  for  (i,j)  = (10, 22)  member  of  the  parameteriza- 
tion set.  Column  e(i)  denotes  expected  values,  while  y(i)  denotes  measured  values. 
The  chosen  parameterization  was  A 6 = 7r / 64  and  A p = h. 


Class 

e{i) 

y(0 

(e(i)-y(i))2/e(i) 

0 

0.37 

0 

0.37 

1 

3.00 

2 

0.33 

2 

12.11 

5 

4.17 

3 

32.29 

26 

1.23 

4 

64.16 

74 

1.51 

5 

101.33 

105 

0.14 

6 

132.42 

162 

6.62 

7 

147.35 

131 

1.81 

8 

142.57 

145 

0.04 

9 

121.76 

107 

1.77 

10 

92.84 

93 

0.00 

11 

63.97 

61 

0.14 

12 

40.12 

40 

0.00 

13 

23.07 

23 

0.00 

14 

12.23 

14 

0.26 

15 

6.01 

7 

0.16 

16 

2.75 

3 

0.02 

> 17 

0.74 

2 

2.14 
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Table  B.2.  Goodness  of  fit  test  data  for  (i,j)  = (55,  7)  member  of  the  parameteriza- 
tion set.  Column  e(i ) denotes  expected  values,  while  y(i)  denotes  measured  values. 
The  chosen  parameterization  was  A 9 = 7r / 64  and  A p = h. 


Class 

e(i) 

y(0 

(e(0  - y(0)7e(0 

< 3 

1.01 

0 

1.01 

4 

2.66 

1 

1.04 

5 

7.05 

10 

1.23 

6 

15.53 

16 

0.01 

7 

29.19 

27 

0.16 

8 

47.82 

39 

1.63 

9 

69.36 

69 

0.00 

10 

90.17 

92 

0.04 

11 

106.15 

102 

0.16 

12 

114.00 

138 

5.03 

13 

112.71 

116 

0.10 

14 

102.90 

96 

0.46 

15 

87.39 

83 

0.22 

16 

69.28 

76 

0.65 

17 

51.48 

46 

0.58 

18 

35.97 

33 

0.26 

19 

23.72 

28 

0.77 

20 

14.79 

9 

2.27 

21 

8.75 

6 

0.86 

22 

4.92 

4 

0.17 

23 

2.63 

7 

7.24 

> 24 

1.20 

2 

0.53 
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Table  B.3.  Goodness  of  fit  test  data  for  (i,j)  = (200,22)  member  of  the  parameteri- 
zation set.  Column  e(i)  denotes  expected  values,  while  y{i)  denotes  measured  values. 
The  parameterization  was  A 6 = 7r / 64  and  A p = h. 


Class 

e(i) 

y(0 

(e(i)-y(i)Y/e(i) 

< 3 

0.73 

i 

0.10 

4 

2.00 

2 

0.00 

5 

5.465 

7 

0.43 

6 

12.42 

6 

3.31 

7 

24.08 

27 

0.35 

8 

40.72 

34 

1.11 

9 

60.96 

75 

3.23 

10 

81.82 

77 

0.28 

11 

99.43 

109 

0.92 

12 

110.3 

96 

1.86 

13 

112.6 

109 

0.11 

14 

106.2 

118 

1.30 

15 

93.18 

79 

2.16 

16 

76.32 

87 

1.50 

17 

58.63 

54 

0.36 

19 

42.32 

42 

0.00 

20 

28.84 

34 

0.92 

21 

18.59 

18 

0.01 

22 

11.37 

11 

0.01 

23 

6.61 

8 

0.29 

24 

3.66 

5 

0.49 

25 

1.94 

1 

0.45 

> 26 

0.86 

0 

0.86 

Table  B.4.  Comparison  of  the  \2  for  the  examples  of  Tables  B.l,  B.2  and  B.3. 


Case 

u 

Xo.95  ,i/ 

x2 

1 

18 

28.87 

20.71 

2 

21 

32.67 

24.42 

3 

23 

35.17 

20.09 
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